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ABSmCT 

Ihis work deals with the dynamical stability of flexible 
cylindrical rods subjected to external parallel ^ uniform, and steady 
flow# A general finite element model is presented for isolated rods 
considering old and new formulations. System matrices are obtained 
using finite element based Galerkin's method and natural boundary 
conditions have been considered explicitly in these matrices along 
with the geometric boundary conditions# The results obtained in the 
form of complex eigenvalues for old and new formulations have been 
compared; wherever possible; with the available results through 
Argand diagrams; and discussed# G-ood agreement is found for pinned- 
pinned and cantilever (with tapered free end) rods except in the ill- 
canditioned range; where even available results can not be relied 
upon# Results for clamped-clamped rods for old formulatioi are 
corrected# New results are shown for clamped-clamped rods for new 
formulation# 

A generalized computer programme is developed to form the 
dynamical matrix; vtoich gives complex eigenvalues by calling a library 
subroutine ZINDIA. An important feature of 1his prcgramme is its 
flexibility. It can accommodate , all the boundary ccoditicns by 
making few changes in the input cards# 
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IHIRODUOTIOH 

1 ,1 imODUCTIOR TO glCW-IEDOOED VIBBATIOH HtCBIEM 

The response of a structure j vshen subjected to an external 
flow and/or containing a flowing fluid, is of considerable practical 
importance and is encountered in various engineering system components, 
such as nuclear fuel bundles (rods) and heat exchanger tubes etc* 

Qlie dynamic behavior of such structures is quite complex and it has 
been a matter of interest for many dynamicists for last two decades* 
This type of problem, where the flow of the fluid, external and/or 
internal, is responsible for the vibration of the structures, is 
termed as "Plow-Induced Vibration Problem"* In reference to the 
orientation of the fluid flow with respect to the axis of the struo-r 
ture , the flow induced vibration problem may be classified as follows s 
1* Parallel flow, i.e*, fluid is flowing parallel to the axis 

of the rod* It is further subdivided as: 

(a) Internal flow, 

(b) External flow, and 

(o) Internal and External flow. 

2. Cross flow, i.e. , fluid is flowing perpendicular to the 


axis of the rod* 


2 


Present work involves the study of dynamic stability of a 
single rod subjected to parallel external flow by the finite element 
method. 

Pluid flowing past a system is a source of energy, which is 
imparted to the system components* This energy is responsible for the 
vibration of the system components which may cause instability in the 
system* This instability, ultimately, may cause the system damage* 

At small flow velocities Instability usually does not occur*' 
As the flow velocity increases, the structure becomes unstable at a 
certain velocity* This velocity is called critical velocity* This 
instability may be of buckling type and/or oscillatory type* These 
instabilities are pronounced enough to cause failure of the system 
components, and hence need to be investigated so that the system can 
be designed to withstand the flow~induoed vibrations* 

Analysis of the flow-induced vibration problems need the 
response, i*e*, the actual motion of the structure, and the stability 
of the structure, i*e«, to know, how sensitive the structure is for 
the dbanges made in the system parameters Involved therein* The sta^ 
bility analysis often supplies sufficient information to considei? the 
vibration problem as solved* Moreover this analysis is comparatively 
simpler* The response analysis requires the complete knowledge of 
the fluid dynamics, i*d., damping mechanism and flow originated 
excitation which may be due to unsteady pressure field or pulsations 


in the flow rate* 
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1,2 HLEVICTJS WORK OH THE DmMICS OB' ROUS STTOJECOIED 
IQ PABAUBl EXTEBHAL Jim 

It was Burgreen et, al |_2 J who first demonstrated the pheno- 
menon of rod vibration caused by the parallel flow of a fluid past a 
cylindrical rod in tbe year 1958, while stud 3 d.ng the heat transfer of 
water flow in g across a bundle of rods in a nuclear reactor core. He 
observed the oscillations to occur at all flow velocities. Also the 


oscillations were found to be relatively constant throughout the whole 
range of velocities which led him to believe that the vibration was 
of self -excited type. He developed an empirical expression for maxi- 
mum amplitude of vibration by dimensional analysis of diffeirential 


equation of motion assuming excitation force to be proportional to the 


1 ' 2 

dynamic pressure p'f ) of the fluid, 
experimental studies p2, 15 , 21, 22^ ^ 


Ihis was followed by several 
with different excitaticxn 


forces, wbsreln the main aim was to study the character of the vibra- 
tions and obtain an empirical expression for maximum amplitude of vibra- 


tion which can be used to solve mechanical design problems in reactors 


and heat excheingers. 


In these experimental studies, various resiparchers assumed 
different forcing functions and obtained empirical expressions for 
macimum amplitude of vibration. Most of them j"3j| selected turbulent- 
boundary layer pressure fluctuation as the forcing function. 

The first analytical and experimental treatment of the prob- 
lem wwas presented by iPaidoussis pi 4-1 in 1966, His mathematical 
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model was based on. Bemoulll-Euler beam theory to describe the flexu- 
ral motions of the cylindrical body) slender body theory for the 
coupled inviscid hydrodynamic farces, and fairly reasonable relation— 
ships for the corresponding viscous forces* He considered three 
boundary conditions; namely pinned-pinned ends, clamped-clamped ends 
and clamped-free ends "With tapered, streamlined downstream end and 
showed that small flow velocities damp the free motion of the system 
but hi^ flow velocities induce buckling and oscillatory hydroelastic 
instabilities* He showed that flutter instabilities are caused by 
lateral frictional forces, and in the absence of hydrodynamic-drag 
effects only buckling instability is possible* Externally applied 
tension was observed to stabilize the system vhile eompresslcn des- 
tabilized it* 

PaidouBsis extended his previous study of the dynamics 

of flexible cylinders in axial flow and presented a general theory to 
account for the small, free, lateral motions of a flexible, slender, 
cyllnderical body with tapered ends, totally submerged in li<iuid and 
towed at steady speed V* Eor particular shapes of the ends and length, 
of tnw^rope, he showed that the body may be subjected to oscillatory 
and non— oscillatory instabilities for V > 0 ; at small V , these 
instabilities correspond to those of a rigid body* He also showed 
that at higher V, the system generally regains stability in the 
above modes » but may be subjected to higher*^ode, flexural oscillatory 
instabilities* ®ae effect of a number of dimensimless parameter cn 
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stability were calculated to obtain the critical conditions of sta- 
bility# It was shown that optimuia stability is achieved with a 
streamlined nose 9 a blunt tail and a short tow-rope* 

In 19739 Paidoussis |j 7 ^ extended the theory by includirg 
material damping 9gravity, pressure forces and fluid damping force, and 
also applied the theory to cylinders in confined flows* In "this work 
frictional forces were taken into account correctly* He showed that 
small flow velocities damp free motions , but at sufficiently hi^ 
velocities hydroelastic instabilities occur. Pinned-pinned cylinders 
are subjected to buckling in their first mode, and then at higher 
flow velocities in their second model at slightly higher flow velo- 
cities ooupledb-mode flutter was observed* Critical velocities fctr 
the cylinder in a closely spaced cluster were much lower than those 
for an isolated one# Cantilever cylinders, on the other hand, were 
subjected to buckling in their first mode and to flutter in their 
higher modes* He carried out the calculations for cantilevers for 
high velocities* These enabled him to observe that the loss of sta- 
bility in the second mode is preceded by regaining of stability in 
the first mode j then in the similar manner second mode regains sta« 
bility at a slightly higher velocity than is necessary for the onset 
of third mode flutter* He also discussed the oase of cylinders in 
axial flow subjected to an arbitrary force field* 

Reontlyln i:; 75 | Paidoussis 20 analysed the stability 


of sULnd 
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Recently in 1975> Paidoussis | 20 J analysed the stability 
of slender cylinders with pinned-pinned ends end cleniped-free ends, 
subjected to axial flow with harmonically perturbed velocity. He 
showed the possibility of instabilities for certain ranges of frequen- 
cies and amplitudes of the perturbations. The parameters which des- 
tabilize the system in steady flow, make the parametric instability 
regions wider. In case of simply-supported cylinders parametric ins- 
tabilities occur over specific ranges of flow velocities. In the 
case of cantilever cylinders the instabilities were found to occur 
in some of the modes of the system and the regions of these instabi- 
lities in Strutt diagrams were quite complex. Moreover, the regioos 
of instability in strutt diagrams became wider as the system approached 
the point of instability in steady flow. 

1,3 ERESEHI. WORK 

In the present work the flow-induced vibration problem, in 
which vibration of the rods is caused by external parallel flow, has 
been solved by the finite element method and the dynamical behavior 
of "the rods with different boundary conditions has been studied. 

Old formulaticn, Paidoussis , and, the new foimulatlcia, Pai- 
doussis » both, have been considered for studying the dynamics 

of rods with pinned-pinned ends, clamped-clamped ends and olamped- 
tapered free ends. The results obtained in the present work have 
been compared, wherever possible, with the available results for 
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old and new fonnulaticaas* Alsoj the effect of various non-dimensional 
parameters on the stability of tho erods has been studied and compared 
with the available results. 

In the second chapter the eq.ua tions of motion are obtained 
in the matrix fozm by the finite element analysis using Oalerkin’s 
method, Eie stiffness* damping and mass matrices for a typical finite 
element are obtained. After assembling these matrices for the vdiole 
domain, natural boundary conditions (as a result of finite element 
based Galerkin’s method) and geometric boundary conditions ace intro- 
duced, finally, the system dynamical matrix is obtained by arrangirg 
these matrices. Method of solution is also discussed in this chapter 

Complex eigenvalues of the dynamical matrix are computed 
using a library subroutine for various boundaory conditions, (The 
various dimensionless parameters, whose effects studied on the sta- 
bility of pinned-pinned rods are : 

(a) effect of dimensionless mass ( 3 ) 

(b) effect of external tension (^) 

(o) effect of internal damping («) 

(d) effect of gravity ( y) 

All the results obtained arc Illustrated and discussed in 
the third ciiapter. lurthenaoire , mechanism of instability is also 
discussed briefly in this chapter. 
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Chapter four contains the conclusioris on the basis of the 
results obtained in the present work. 

Further, in Appendix A, various fluid forces resulting from 
the fluid flow have been discussed iu detail'. Also governing eq.ua- 
tion for small lateral motions for an isolated rod has been derived. 

Appendix B deals with the old formulation of the problem 
and its finite element analysis. The difference between old and new 
formulations has been discussed in brief. Appendix C gives the boun- 
dary conditions for a cantilever rod with tapered free end, 

lastly, the computer programme has been introduced in 
Appendix D whose detailed listing is given in Appendix E* 



OHAEEBE II 


PaRMJIAIIOH Qg ERCBIEM AID OOMPUIAglOKAli APIROAce 


In the present chapter the fifth order differential equa- 
tion govemiog the transverse moticjn of a rod, in an infinite fluid 
flowing at a velocity V is converted into simultaneous differential 
equatio 2 as by finite element method using G-alerkin's method* Purthery 
these equations are rearranged to fom standard eigenvalue problem 
to obtain complex eigenvalues* 


2,1 FINIIB BIEMBHT MALYSIS ; 


!I!he governing equation of motion for the small transverse 
vibration of a flexible slender rod of circular cross-section immersed 
in an infinite fluid flowing with a uniform velocity V can be 
written asj Paidoussis HtI j 


VI 




3x^ 9t 


34 

+ BI — ^ + 


3^ 
G m„ — ^ 


+ 2 


0 mV 
m 


9 y 

■f ' 33 c 3t 




2 3 


-JL - { 5 
9x 


[? + (l-2v ) (5.)] + 


+ 0^ + t (1- “Is ) - X ] 



Pf 3>V 0^(1 -«) Pf 

'’f ® “b It * “r 


s= 0 , (2,1) 
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with the boundary conditions as 


5 4 

(t + (1 - 2v ) ( pA )) ^ 

3 x at 


+ K^q y=0 ory=0 at x=:0, 


BI + p I K |-J=0 or~^=0 at x=0, 

3x ax at 

El + pi -^-=2L-„ -.(!!! + ( i >. 2 v ) (pA 

ax^ ax^at 


- y=0 or y=0 at x=Ii, 


El ^ + U I - - ' I — + K = 0 or = 0 at x = L, 

ax at 


ax 


ax 


( 2 , 2 ) 


inhere Z and W are given by Bqns. (2,3) and (2,4), respectively, 
and 


p sr internal damping coefficient representing a Kelvin-Voigt 

type of damping in the material of the rod» 

BI a flexural rigidity of the rod, 

0 e coefficient of virtual mass, 
m 

m^ = virtual mass of rod per unit length. 

m ss mass of rod per unit length, 

r 

V ss mean flow velocity of the fluid. 

6 = equal to 0 or 1 accordingly as the rod is free or not 


to slide axially at x = I# 
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I S3 externally applied, uniform tension. 

•yj S3 Poisson’s ratio, 

p = steady-state pressure of iiie fluid. 

A ss cross-sectional area of rod* 
pA ss pA at X = . 

S3 density of the fluid. 

I) ss diameter of rod. 

0^ « skin friction drag coefficient. 

Ojj 53 coefficient of viscous damping. 

0^ 53 base drag coefficient. 

It *5 length of rod. 

X SS axial co-ordinate, 
t 33 time. 

y S3 lateral dtsplacement of rod. 
dO i 

yf. SS displacement spring constants at x sr 0 and x ss L, respectively, 
^dli 

^tO 

*33 rotational spring constants at x ss 0 and x = It, respectively. 

^tli 

Por the sake of simplicity, various terms in Eqn. (2.1 ) 

2 

associated with ^ and are put in the following manner : 

Z =3 (C^ -{6 [ T + (1 - 2 v) ( ^ ) ].+ [-1 Pf 3) 0^+ (m^-3a^)gl 

1(1 --^6) It] +[“ Pj 0^ (1 - «)]}), (2.3) 
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(m + C m ) + y I — - +20 in„ V 

Tf'TriT ^C ifl. ryi-PfivC^ 


m I ^^2 ‘ . 4 . 

3t 3x 3t 


m f 3x 3t 


+ X 


8>5 __3V^> . 3M^> ... .3 


3 t + El + z — 5=J- + I (fe^+x) -Sj— + ff- ) 


3 X 


3x“ P 3x 

1. = 1 y^9 • ■ • • • • C2*S^ 

Now, 1316 above equation is integrated by pai'ts to introduce 
the natural boundary conditions and reduce the order of derivatives 
simultaneously. 


dxj=0. 


Second term is integrated twice and it 3 results, 


li ^(e) 

/ N UI-~^ 

0 3x^ 9t 


dx = 




4 (e) 
Z_ 


Qx' 9t 


h 


dN. 

1 

dx 


2 

h d N. 

^ yi^ 

0 dx 3 X 3 1 


33.(=) 


dx 


yi 


3^M 


3x^ 3 t 


h 


(2.9) 


Third term, after integrated once, yields, 


f N 2 0 m V 
J - m f ^ 3x 3 1 


2 (e) 

TT * ' 


dx = 


u: 


(e) 


N. 2 C m » V _ 
X m f 3t 


-/ 


h 

0 


H aN (e) 

r-— 2 0 m V dx . 

Q dx m f 3t 


(2.10) 


/ H El ^ 
0 ^ 3x 


Pifth term is also integrated twice and it gives, 

34 (e) I 35„(e) id 1 32,,(®) I ^ 


dx = 


N^EI 


JL 


3x" 


dN. 

3x^ 


ii a^H. 


- 32 M 

-J-BiiJtg-i 
0 dx 3x 


dx . 


(2,11) 
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Sixth teim reduces to, 

h 2 ^ /V 

r HT rr 3 ^ I h dSr 

j U Z -—I ax = IT. Z _J, 7 3,/ , 

0 8x^ 1^3^ 0 0^ 3^^ * 

Seventh texe. tume out to he''^ totSiet^eSlt”'"^ ^ 
integration by parts gives, 

/ N. W((x + x) -^5^4. 3.7^®\ . K , ^ av^®^ 

0 ^ P 9x2 3^^ dx = [ IT. w (xp H- x) |X-~ 


h dU. (e) 


0 


0 

(2.13) 


Substituting Eqns. (2.9), (2.I0), (2.II), (2.12) and (2.13) in Eqn. 
(2.8), it takes the following fom. 


W eUl y JLV. 

m f 13^2 

3„(e) d*u. ,2 (e) 

^ \ 9 t ^ ^ 


2 a 2 

dx ox 


p n V XU V * "*' ifci- 

te 3x2 9 1 m f dx at 


. (e) 


dx 9x 


<pr. 

« 2 -1^ 
dx 

3x 

- Wx 

] 

9 4 (e ) 

miiim 4^ 

h 

dIT. 
1 

3x^ 9t 

|o 1 

dx ' 


■p cbc ax 

3 (®) .h 


y I 


ajL 

3 x 2 3 t |o 


3v'®> " |«l,8V=) |h 


\ 2 V + IT El 

10 "" ax^ 


+ H. Z 
X 9x 


(e) 


0 


h 


0 


+ I W (xp + x)|^- 


(e) 


dx - 2 

a X 


* 0 


= 0 


i — 1 , 2 ,«ee*, f 


(2.14) 
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Agato, to have -Hoe ayatem mtrioes , Eqn. (2,7) le substituted 

into the above Ega. ( 2 .I 4 ) sad aearraagiag, one gets for the sleaeat, 
e, 


[m] o] { 4 }“® 




If. yl ^ 
* a ^ 


94 (e) jh 


Ox t ' 0 

4„(e) ,h 


w m 
^ 3 x^ 





m Y 

^ lo 




(e) |h 


9y^®5 1^ 

JH Y ^-4 

^ lo 


t a^^yCe) h 
El 

9 x 0 

, g 2 (e) , i 
N, El 2-5 


, gaje) I h 

“f .- 1 — „ 

3x ' 0 
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H, 


1 3x 

3Z 


Z 


h 
0 

(e) , h 


2 “ 3x 




Z 




(e)| 


f 3x 


0 

h 

0 




V*-. 

where 




tml 


(e) 


r 


K. W (x + x)^ 
1 p ''3x 


(e) 


1„ W (x + x) 

2 ' p 3x 


h 
0 

(e) , h 






iT„ W (x + x) 
f p 3x 


(e) 




= mass matrix for the element 


0 


h 


>= fo} 


( 2 . 15 ) 


/ (m^+O^mP {H} L'S'J ite 


U) 


[c] ^ total damping matrix: for the element 


[Cl]^®^ +[c23^®^ 


[ c3l 


(e) 


,(e) 


[ cl ] = internal damping matrix for the element 


h 


n If 


= f u I {If } L ^ J ^ 
0 


i(e) 


[c2l =s coriolic force dampirg matrix for the element 


h 




( 2 . 16 ) 


( 2 . 17 ) 


( 2 . 18 ) 


( 2 . 19 ) 


f 1 Ce / 

I c3 J = viscous drag damping matrix for the element 


= / X{]iI}|_]jJdx 
0 


( 2 . 20 ) 
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[k] 


(e) 


s= total stiffness matrix for the element 


= t - Ck3]^®^ 


[k1] 


(e) 


s= elastic stiffness matrix for the element 


( 2 . 21 ) 


=/ Ei{ I } LnJ ax 

0 


( 2 . 22 ) 


(e) 


[ k2] ^ = stiffness matrix for the element due to curvature term 

Lz 


9x 


= / (2 + Wx){H}Lirj dx 


(2.23) 


0 


and 


,(e) 


[k3j ^ = stifjEness matrix for the element due to variable 
cuivature term x 


3x 

k I t 

= / M sii IS }LJJJ<3x 
0 


(2.24-) 


.(e) 


To calculate these matrices one has to choose 7' '. The 

(e) 

following canditions must be fulfilled while selecting this y 
to eaasure monotonic convergence ^ Huebner I 9 1 * 

(i) Completeness j i.e.j all uniform states of y ^nd its partial 

derivatives upto the highest order appearing in Bqn. (2.14) 

(e) 

should have a representation in y \Ah.en in the limitj "the 
siae of the element shrinks to zero. 
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(ii) Ccmpatibility , i.e*, at the inteirfaces of the elements y 
and any of its partial derivati-res one order less than the 
highest order appearing in Sqn® (2.I4) must be continuous. 

The following shape functiouj usually used for beam problems, 
in terms of generalizsd co-ordinates has been chosen for the present 
problem, since it satisfies the above requirements of completeness 
and compatibility, respectively, 

LirJ = L5 i^(3-2£i) Ejt e / (3 - 2 5 p --t, e/ bj , 

.... (2.25) 

where = 1 - " and ^2 " h * 


Introducing the shape function into the Eqns. (2.16), (2.18), 
(2.19), (2.20), (2,22), (2.23) and (2.24) and integrating, one obtains 
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Ic2] 


/V 2 C Y 
V® ) in f 


60 


[c3] 


(e) 


Xli 

420 


, ^ (e) 

[W] 

h 


[k2] 


and 


/ ^ (Z + W X ) 

Ce) -P-1 


3011 


30 

-6h 

-30 


6h 30 

0 6h 

-6h -30 


-6h 

-h" 

61i 


36 


3!a 

411^ 


-36 

-3li 

36 


[ k5l W_ 

= 30 

I 

Symmetric 

DEherefore Eqn. (2.15) 1)60011165, 


Symmetric 


18 31i 


-18 

-311 

18 


6b 

L- 


-6b 

0 



""l56 

22b 

54 

-13b 


CM 

13b 

-3b^ 



156 

-22b 

Symmetric 


CM 

12 

6b. ' 

-12 

6b 


2 

4b 

-6b 


, 


12 

-6h 

Symmetric 


CM 


31i 

-b^ 

-3b 

4b^ 


0 

-bVc 


0 


3h 


(2.28) 


(2.29) 


(2.30) 


(2.31) 


( 2 . 52 ) 
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Sot the #io3e domain of the problem, Eqn. (2.33) becomes 
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- 2 C m^ V 


is: 


(e) 




I 

< 


xa t at 

0 
0 




(e) 






+ Z 


0 
0 

• 

0 

3 „(e) 


> 




0 

3 y(® ) 

2 ^a”t~ 


0 




0 

0 

0 




JC 


Jmrf-I 


V' 


m li 


lx 


(e) 


(ei 




VI ^’.)r 

Sx"^ 3t 


'w 


m+l 




0 




(- EI + 111 -^) 


3x 


3x^3 t 


g3 (e ) j. ^ 


9x" 


3x^ 3t 




0 


m+1 




(ei 


2 (e) 


ijSL 




3x" 


0 

0 

+ U I 


2 . 

3x 3 t 


m+1 


(2.34) 
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wfcere 

[ Ml— = assembled mass matrix for ttie vdiole domain before applying 
the boundary conditions 

[ C]_ = assembled damping matrix for the whole domain before apply- 

-DD 

ing the boundary conditions 
and 

[ = assembled stiffness matrix for the whole domain before 

applying the boundary conditions « 

For Gonputational ease, Eqn. (2.34) can be rearranged as 
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.\ie) .4 (e) 

- (El ^ ^ I 


8 X 


8z' 




0 

0 






0 


-(ei 


82 (e) 


+ y I 


|3je) 


dx 


9 x at 




(BI 

9x-^ 


3 (e) 

TT ' * 


0 


+ y I 


(e) 

8_|_) 

9x^ 9t 


m+l 


(ei 




0 

0 

* 

0 


8 Me) 


i 

y 


+ yl 


.3 (e) 
9 y 


2 ^ 
9x t 


mi-1 


(2.35) 


In Eqn« (2.35), the values - 2 0 m _ TT and + 2 C m „ Y aa?e 

la X mi 

added up with the elements 0^ ^ and 0„ - - ^ of the assembled 

I 5 I I ^ ♦ 

damping matrix, respectively. It should also be noted 14iat the co- 
efficients of 9 and 0 are simmed-up with the elements 2 
^aiH-1 2m+2 assembled stiffness matrix, respectively. 

if ter substituting the boundary conditions, Eqn. (2.35) can 
be written as 

[M]{ V}“ + [0]{y}“ + [Kj {y}"^ ={o}, (2.36) 

where ( M ] » t 0 ] and [K] are the assembled mass, damping and 

stiffness matrices, respectively, for the tdaole domain of 1iie problem, 
after applying -the boundary conditions. 

Equation of motion (2.1) derived by Paidoussis t17l is much 
moire general than that derived earlier, Paidoussis [14] 


t 


as it 
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takes iato consideration liie gravity and pressurization effects. 

In addition, it eliminated the error of nan-inclusion of the fric- 
tional term ( 9y/8x) in the earlier woik. Earlier formulation 
[ 14 ] is refer3?ed as "Old Pomulaticaa" and the later f ouQulation 
t 17 1 used in this chapter as "Hew Pormulation" in the present 
work and are given in Appendices A and B. 

The author found it fascinating to solve the old fonaulation 
problem also by the finite element method and is given in Appendix B. 

2,2 GOIfflEUTATIOHAl AEERQACH •. 

The solution of Ihese hanogeneous set of differential equa- 
ti033S is obtained as described in Ifeirovitch fill • As per th-ia 
reference the inertia, stiffness and damping matrices should nece- 
ssarily be symmetric to arrive at the solution. But it is not essen- 
tial for the matrices [ M ] , [ K ] and [ C ] to be symmetric , refer 
Praaer, Duncan and Collar [6 ] . 

Equation of motion of forced vibration of a damped system 
is 

CM 3 { y (t) } + [ 0 ] { y (t) } + [ K 3 { y(t) } = {Q(t)} , 

.... (2.37) 

Equations (2.37) repiresent a set of coupled ordinary linear 
diffei«ntial equations vshich is difficult to decouple unless matrix 
[ C ] is a linear combination of matrices [ M ] and [K] , vhich can 
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occur in particular cases only, lor a general case, equations (2.37) 
can be 37educed to a system of first order diffearential equaticaas as 
follows , Meirovitch [11] , 


Premultiplying Eqn. (2.3?) by [ M ] 


-1 


ry (t) > +[ M 1“^ [C] { y(t)} + [M] [K]f y(t) } = E M ]“^ {Q(t) } 


-1 


,-1 


let { q(t)} =. 


(K-b)} 




and 


*> 


[ 2n X 1 ] 


r{o 


( 2 . 38 ) 


J 


][ 2n X 1 1 

(2.59) 

This leads to a set of 2n first order ordinary differential equa- 
tions 

tire q(t) } + C K 1 {q(t) } = {y(t) } , (2.40) 

where 

[To] [mTI ^ -*[ M] [ 0] 

and [ K ] = (2*41) 

_[ 03 [ Kl^ 

are real matrices of order 2n. 

Prom Eqn. (2.40) we obtain the homegeneous set of differen- 
tial equations in matrix form as 

til { q(t)} + fKJ{q(t) > s { 0 } (2.42) 

The solution of Eqn. (2*42) is obtained by letting 

{q(t)} =. e^'*' { } , (2.43) 

Where 0 ^- represents a vector consisting of 2n constant elements'i 
Eqn. (2*43) j when introduced in Eqn. (2.42), leads to the eigenvalue 


problem 
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[i ] {^} + [ ic ] {jzf } = { 0} , 

which can be rewrittien in liie form 
[D] { jZJ} = ~{ jZi} , 

A 

where 


D = - 


[ K ] [ M ] 


CM] 


-1 


Col 


[0] 


Co] 


- [Kf 


Co] Cm] 


Cm] [ c 3 


C i3 


L [K]"^ Cm! - CkF^ Cc] 


[2n X 2n] 


(2.44) 

(2.45) 


(2.46) 


plays the role of dynaaical matrix. The multiplication o£ parti- 
tioned matrices is handled as if the submatrices wei« OTdinary ele- 
ments. The identity matrix [ Ij in Eqn. (2.46) is of order n. 


The eigenvalues of the above dynamical matrix (2.46) will 
in general be complex. Real parts determine the stabilily of the 
system and Uie imaginary parts give the frequency of vibration. If 
the real part is negative-, ttie system will be stable and vice-versa. 
The system will have oscillatory instability (flutter) if real part 
is positive and imaginary part is non-zero | whereas if the imaginary 
part is zero, the system will undergo buckling instab iliiy. 

rp ba formulation presented in this chapter is general and 
proves to be versatile because the effect of various boundary 



conditions (geometrical and natural) can be studied easily. Also, 
the effect of various system parameters such as velocity, axial ten- 
sion, fluid pressure, gravity etc. on the dynamics of rods can be 
studied with special interest* 

A computer programme has been prepared to compute the eigen- 
vaHues using built-in subroutine ZIHDIA* lasting of camputer programme 
is given in Appendix E* Computer programme has been discussed in 
Appendix 



CHA.P!EER III 


EESUIffS MD Discussion 


'■ The results obtaiiaed in the form, of dimensionless parameters 
are included and discussed in this chapter. Various boundary condi- 
tions of the rod and the effect of rariotts system parameters on the 
stability of the rod are illustrated in detail. Complex frequencies 
of the lowest four modes have been obtained as a function of the 
dimensicnless flow velocity and Argand diagrams have been plotted. 
All the results are given in temns of the dimensionless parameters. 
Lastly, the mechanism of instability has been discussed in brief. 


50 DliiEjaSIOUIESS PARAI^EITERS 

following dimensionless parameters have been used in this 


work, 

3 = 

u = 

r s: 
a = 

Y = 

T) 


Dimensionless mass = m^/(m^ + 12^,) 

Dimensionless flow velocity = (m^/Bl) VL 

— 2 

DlmensionlQss axial tension = T L /El 


Dimensionless internal damping coefficient 
= {1/ [ E (m^ + m^) ] } V / 


Dimensionless gravity parameter 
Dimensionless pressure parameter 


_ (m^ - m^) g T?/B 1 
= { ) L^/EI 
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1 /2 

Ee(oj ) s= Real part of the dimensicaaless frequency ={ (m^ + 


Im( to) =: Imaginary part of the dime33sionless frequency 

1/2 _2 

** _ -t- Hi / ri. 1 r 

f r 


a 


a. 


"f 

°b 


= { (m„ + tOxlj'' 


3 

s= Dimensionless displacement spring factor = K^D /BI 

2 

s Dimensionless rotational spring factor = K^D /BI 
= Dimensionless viscous drag coefficient = CjjL(m^^/Bl) 
s= Rriction drag coefficient = (4/ tt ) 0^ 
s= Base drag coefficient = (4/'Jf ) 0^ 


1/2 


3 • 2 CQWER&ERCB OF GHB EIGBFyAIDES : 

Finite element analysis says that as the number of complete 
and compatible finite elements is increased^ "tiie results converge 
monotonioally to the exadt values. However, one has to choose finite 
number of finite elements. For the present work, first for every 
boundary conditian results were obtained for a few ncn-dimensional 
velocities for various numbers of the finite elements. It was obser- 
ved that for five elements results were quite acceptable for all the 
cases. Sample results for pinned-pinned rod and cantilever rod are 
shown in Table 1 and 2,. for WI^ frcm .2, to 7 m&.frm. Z to 5,, febpec- 
ti've3y,(pp,30,5l)jtt should be noted that ttie change in the frequencies 
after NED s= 4 is negligible. Therefore, it was decided to use five 
finite elements for determining the results. libreover, this saved 


computational time also* 
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3.3 OOMPIEX FREQUENCIES OP ROIS ; 

The dimensionless complex frequencies obtained for the pre- 
sent problem are plotted in Argand diagrams for different boundaiy 
conditians of the rod* The imaginary part of the dimensionless 
value Im(a) ), which represents the frequency cxf vibration is plotted 
on the abscissa and the real part Ee(ti) )5 on the ordinate* The various 
parameters involved in the problem are ot, ti, u, ® t ’ 

°f » % °d- 

3.3*1 PIMED-PIIMBD ROD ; 

The boundary ccnditicns for pinned-pinned rod are 


y = 0 at X = 0 and 

X — 1, 

(3.1) 

El ^ + U I = 

s 0 at X ss Oj 

(3.2) 

3x 3x 3t 



El ^ + yl 

=s 0 at X = 1. 

(3.3) 

3x 3 X 3 t 




Since pinned-pinned rods are of considerable practical inte- 
rest as compared to the cantilevered rods and clamped -clamped rodS| 
this boundary condition is studied in detail and the effect of various 
dimensionless parameters on the stabili'ty of the rod is discussed* 

3*3. 1.1 OH) FORMJIAT I QH: 

In this old formulation, transverse ccmponent of tangential 
fractional force is neglected, see Paidoussis [ 14 1 and Appen- 
dix B. Pig. 3.1 shows the dimensionless complex frequencies of the 
lowest three modes of a system as a function of the dimensionless 
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flow -velocity u , for 6 = 0o1,ec^= 1, F = 0, “ These results 

•were compared -with -the results of Paidoussis [14 ] and it was ol^served 
that results were matching v/ith those of available resul-ts except fcsr 
certain velocity ranges, i.e^, in the first mode for u ^3 and and 
in second mode for 5< u< 6, Careful study of the available results 
indicated that there was ill-conditioning for the above mentioned 
velocity ranges, i.e», for small changes in u, a large change in 
complex frequencies was obseived. This led the author to believe 
that there is ill-conditioning in the system in these velocity ranges 
and the convergence is incomplete for five elements, therefore, the 
results for these velocity rarges, -vdiere ill-conditioning was believed 
to be present, were studied by increasirg the number of elements from 
2 to 8'« These results are tabulated in Table 3» 

Still for 8 elements there is seme difference in the present 
and Paidoussis’s[ 14 3 results in these velocity ranges. Moreover, 
the results of Paidoussis [ 143 may also not he reliable for the 
ill-conditioned ranges, 

3 . 3 , 1 . 2 lEW PCEMUIATIQlSf i 

Case (a) : 3 = 0,1,ec^= 1,6 = (^=: 1, = Y“ ri == = 0, 

Pig, 3,2 shows the dimensionless complex frequencies of the 
lowest three modes of the system as a function of the dimensionless 
flow velocity u. The results match very -well with these ob-talned by 
Paidoussis fl7 3 except for certain ranges of u, where ill-condition- 
ing was observed to be present, i,e,, 6 <. u<, 6,5 second mode 



34 



LTN 

m 

LO 

(1 

/ 

V 

di 

I 

ty V 

4 

j 

ljr\ ON cj C>- on r- 3 r\ 

0 C\J v£) vO ON vO J- 

ON ON c\i 0 ON On On 

• • • • • • • 

T- 0 0 0 CN On 0 

OJ (M (M CM V- •«“ T- 

3 

w' 

0 

3:1 

m r- CJN CM cn cn 

^ 0 0 t- T- •r- •?= 

-r- 'T- ^ *r“ T-* 

• * ♦ • • • ♦ 

0 0 0 0 0 0 0 

p{ o 

0 | r- 

M • 

5 o a 

.1 s ■ 

^ ca ^ 

. “T 

ca ^ 5 

^ § p 

1 S p 

^ p ^ 

a 1 § 

6 s ^ 

a 9 

a 1 

^ H 

S Ah I 

! 

1 

ijr> 

J" 

• 

f! 

PJ 

Q 

3 

trx 

CO 

J- 0 0 0 0 0 0 

• 

0 

C 

d 

4 J? 

U] 

H 

"a 

0 

p:; 

0 CO Q\ v- ON onoN oc^ ono J“t7 

S' Slfti, 55 83 ! -iS -SI: 

° °? "T? "79 "79 " 7 ? " 7 ? 

T-* 

# 

on 

" 

?5 

0 

3 

>w^ 

B 

H 

\r\ 

on 

CNO 0 0 o"0 0 

* 

0 

w 

0 

H 

3 

0 

pel 

0 O^D vDlO\ OnOO y-On Q>nO 

ON 0 ON 4 - + OnCO 0 \^ 00 ON 

CM 0 tr\ O-T** Cn-y— EN-^ IN-y* Cn-y— 

• •* •• #• •* •• •• 

99 . 99 99 99 99 9? 1 

CO 

on 

« 

on 

n 

0 

■B 

3 

1 

g: 

S V- 0 0 000 0 

• • 

0 0 

s 

•p 

0 

H 

'3 

0 

fci 

0 <X> CM 1 >- QO t>- ITnCO cm r- 

ON CJN CM Cn. OJ CM J-.J* 'iXNljCN tTNUN 

CM CM Otrs VOO ^0 vOO M 30 

• • •• •• •• •« •• 

0 0 00 OQ OQ OQ 00 

It I + I ^ t + 11 

0 
cr 

« 

O'" 

n 

r^i 

1 r*- 

0 

B 

s 

P 

cc 

! H 

3 

s 

V" \CN 

■UN 0 

CO m 0 0 0 0 0 

• • 

0 0 

i-o 

Ph 

CJN ad CMCN ^\D ONf- ^00 MD^ 

CO CJN o-Cvl So.^ vr\j- \iMrN J-i?n 

CM CM r-J- OVN 0^ OtfN OVN 

0 C* «9 •# •* •• 

0 0 00 00 00 00 00 

1 .1 1 t t 1 ! I 1 1 

CZi 

or 

T- 

i 

or 

11 

r- 

r" 

■ g; 
. ^ 
1 c 

“ g 

4: 
u 
5 h 

!c 

' \0 CM H- ON Xtn 

j CO CJ\ T- -r- 1 £N 

CO cn X- 0 0 0 

0 C 0 

00 0 0-0 

5 

3 y- 
^ •= 

a 

pc 

A ns CO CJN 0 0 OCO OTcM 

s ro CJN ON 0 0 ^-4^ 5 :?':^ 

- CM CM CM rO , rO. cM CM cn 

• p « 

5 0 0 0 0 0 C 3 CD CDCD 

§ 

rC 

E 

r 

' 

0 

-I" 

§ 

5 

3 

-1 

3 CM CO jd" ICn sO OO 

Ji 

1 . ^ — 



55 


and 7 i 10 in. the third mode, Studyiig Pig, 3*2 and Table 5 it 
vras thought that not much t/ill be gained by increasing the number of 
finite elements* 

It is observed that small flow velocities damp free motions 
of the rod* It is also noted that the behavior of the pinned-pinned 
rod is changed a lot by the corrections made in the frictional force 
terms* In the older fonnulaticn, Fig* 3*1 » buckling occurs in the 
first mode only and flutter occurs in second and third modes* ’ithereaS} 
in the new formulation, Fig* 3*2, buckling can occur in aU the three 
modes and there is coupled mode flutter at uas 6,48 (first mode for 
u > 6*48* It may be noted that there is no coupled-mode flutter for 
u > 6*48, as wrongly interpreted by Paidoussis [ 17 ] • 

Oase (b) ; 8 = 0.48, ec^ = 0.25, 5= = 1 , c^ =a= y= n=r = 0. 

Fig. 5.3 shows the dimensionless complex frequencies of the 
lowest three modes for the system as a function of the dimensionless 
flow velocity u* Considerably good agreement of the results with 
those obtained by Paidoussis [ 171 is observed* In the first mode 
ill— conditioning is present for 3^ u<3*25 and 6*25< u< 6*4* Slight 
discrepancy is observed in the results for 6,25 <u <6.4* 

Results for 8= 0.1 and keeping the other parameters same 
are also shown in Fig. 3.3 for u upto 3.1* It is interesting to 
note that the increase in 0 tends to increase the damping in the 
system for small values of u. In fact, it can be seen from the 
coefficient of from Eqn* (2.6). 
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In this case second mode does not cross Im( O) ) - axis. 

Thus no buckling occurs in the second mode. The first mode branch 
first crosses Im( to) - axis through the origin at u k 3.H and again 
at u s! 6.29* It leaves the real axis at ua 6.295 and flutter occurs 
for u ^6.295* There is no coupled-mode flutter, as iivrongly inteiv 
preted by Paidoussis [171 * 

Case (q) ; EPEBOT OP TENSION ; 0 s 0.1 » e = 0.25, T = 10, 

«^=:a=Y = n==0. 

Pig. 3.4 shows the dimensionless complex: frequencies of the 
lowest three modes for the system as a function of dimensionless flow 
velocity u. Par this case, Paidoussis [I 7 I obtained the values of 
first buckling velocity, second buckling velocity and the critical 
velocity far coupled-mode flutter as u^^ = 4.455, u ^2 = 7.04 and 
ss 7.11. The author obtained these values as u^^a 4.456, 
and 7.12, This sli^t discrepency between the results is due to 

ill— conditioning in the system near thesQ velocities* It is observed 
that extejrnally applied axial tension has stabilizing effect, i.e,, it 
raises the critical flow velocities. 

Case (d) ; EPPEGT OP IHTBBEAL DMPIEa : 0.1, c = 0.25, 

o = 0.003, «= 0^ = ■'5 =r= 0. 

Pig¥3.$ ^ows the dimensionless complex frequencies for the 
system as a function of the dimensionless flow velocity u. TIm 
dynamic behavior is very similar to the Case (A), see Pig. 3*2. 



37 


The first and second buckling velocities s as obtained by Paidoussis 
[173 are = 3.143 and 

for coupled-mode flutter is u ^ = 6.52# In the present case, the 

cO 

author obtained these velocities as - 3.141 

6.35 • Slight discrepency in these results is due to ill-condi- 
tioning near these velocities. It is noted from Pig. 3.5 that coupled- 
mode flutter occurs at u a 6.35 and flutter occurs in the first 
mode for u > 6.35. 


%2 


~ 6.29 and 


u ^2 s: 6.28 and the critical velocity 


Case (e) : EPPEOT OF CRAYITY : 0 =0.1, e = 0.25, y= 10, 

= °m “ =« = Ti = r = O. 

Pig. 3.6 shoi/vs the dimensionlesa corapler. frequencies for the 
system as the function of the dimensionless flav velocity u. The 
first and seccaid buckling 'velocities obtained by Paidoussis [17 3 
are Uj^^ = 3.1Q7, u^g = 6.30 and the critical velocity for coupled- 
mode flutter is u ^ = 6.44. The author obtained these values as 
“3.2, u ^2 %0 “ Here also ill-conditioning 

is present for the velocity rasages 3*1 < u < 3*5 and 6 < u <6.5. 
It is observed that the inclusion of gravity parameter subtly 
reduces the values of Im(co ) while the effect on Re(u) is negli- 


gible • 
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3*3.2 G IAMESD-CLilMEEP ROD : 

The boundary conditian for clamped-clamped rod are 
y = 0 at X c 0 acid x = I, (5*4) 

y’= 0 at X = 0 and x = 1, (3*4) 

3'.3*2.1 OID PORMJMTIOIjf : 

Pig* 3*7 shews the dimensicoalfiss cemplex frequencies of the 
lowest three modes for the system as a function of the d i m en sionless 
flow velocity Uj fesr 3=s0.1,ecj=l3^ =0'* Results 

obtained by Raidoussie [ 14] are also shown in the sane figure* It 
is observed that the results do not match. In this case no ill- 
conditioning is observed in the present work or Paidoussis's work 
[I 4 I * Author’s results are correct as the deck ■ used for pinned- 
pinned rodj modified for the clamped-clamped case was used* Correct 
results for u = 0 confirmed that the boundary conditions were 
applied correctly. 

It is observed that buckling will occur at u“ 6*25 iJi ‘tk.e 
first mode which contradicts the earlier result, Paidoussis [14] • 

Buckling does not occur in second and third modes. Plutter instabi- 
lity occurs in second mode at u “ 7*5 and in third mode at u« 10*7* 

It is interesting to note that third mode regiens stabilii^r at uat 11*6* 

3 *3*2. 2 mm PORMJIATIOR ; 

Pig. 3.8 shows the dimensionless emplex frequeijcies of the ■ 
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lowest three modes for the system as a functim of the dimensionless 
flow velocity u, for B = O.IjE = 1,6 = 0^ = 1 , =a=:Ys=ri==r 

= 0* Ccmparirig this with Pig. 3«7 it is seen that due to inclusion 
of transverse component of frictional force, 5*^ , the behavior 
changes considerably. low instability occurs by buckling in its 
first mode at u « 6,25 and in the third mode at u « 12, And it is 
interestirg to note that the3re is no flutter instability when new 
formulation is used. 

3.3.3 OAmilETBS ROD WITH OlAIEBED 3i!REE ETO s 

The boundary conditions for a cantilever rod with tapered 
frea end are, see Appendix 0, 

y=y’ = 0 at x = 0 (3.6) 

BI -^ + U I + f m . V ( ^ + 7 1^) 

3x^ 3x^9t ^ ^ 

3 2 

- 0a„ + f m ) x'g ^ = 0 (3.7) 

r I 9^ 

32 33 

aaxi BI • + V* X — ? O, at x Ii , (3.8) 

3x 9x 3t 

1 

vhere x% » r* J A(x) dx. 

^ 1^1 

3. 3. 3.1 OID BaRMJIATIOIir i 

Pig. 3.9 shows the ddmensiculess complex frequencies of Ihe 
lowest four modes for the system as a functions of dimensionless flow 
velocily u, for g= 0.5, 5 = ^ ® 0.01, 



r=: 0* Eesults match completely with those obtained by classical 
methods, Paidoussis [14 ] » Both the results are shown in Pig» 3»9 
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{ 

and'f 


Here , it is important to note that the column matrices 

I I . .,(®) 


Iff 2 m V 


33 ^ 


0 






N. ( El „ 


4- y I 


4,, (®). 




^i 


JL 


3 

gx^ t 


} found in Eqn. (2* 14) must 


XsL 


be taken into ccnsideraticn in mass, damping and stiffness matrices, 
as shown by Deb [ 26 3 . In the analysis these can not be neglected 
as suggested by Szabo and lee [24] * 


3.3.3.2 MEV FQRMnAIIOl : 

Pig. 3.10 shows the dimensionless complex frequencies of the 

lowest four modes for •tiie system as a function of dimensionless flow 

velocity u, for 3= 0.5, e ” 1, 6s= 0, - 0.8, = 0*01, 

c =a:i)ri=: Y=s c =0. Complete agreement was observed with tiie s^esults 
b d 

of Paidoussis [173 upto u = 8. After u> 8, ill-ccsnditioning was 
observed and for this velocity rarge, the convergence was studied by 
increasing the number of finite elements from 2 to 9. Eesults are 
tabulated in Table 4. Also in Pig. 3.10, for u > 8, results with 
HELss 5, EEIi s: 8 and the results obtained by Paidoussis [ 173 are 
shown. Satisfactory ccnvergenoe is not observed for u> 8 because 
of ill*«conditicining. 

Comparing the Pigs. 3.9 and 3.10, it is observed that there 
is no basic changs in the behavior of the system by including the 




-^.692 16.05lf 
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transverse component of the tangential fricticsnal force in the 

new fonaulatian* Buckling occurs in fiist mode at u^a 2.04 and 
flutter occurs at u^^ « 5*15 in the second mode and at '^^ 2 “ 
in the thiii mode. 

It should be noted that the results can not be relied upon 
for velocities u > 8, as ill-conditioning occurs in this velocity 
range. 


3.4 'PHB MBCHABISM Qg IHSOifeBIIlIY : 


After studying the various results f the mechanism of various 
hydroelastic instabilities for rods subjected to external flows may 
be explained in brief as follows*. 


5 , 4,1 BUCKIMO IBSIABIIITIES ; 

2 8^v 9^V 

It can be seen that the term m_ 7 — ^ in 2 ^ of the 

^ 9x: 9x 

governing equation of motion, Eqn« (2,6), behaves more or less as an 
equivalent compressive force ishioh increases with the increase in 


flow velocity I thus reducing natural frequencies of the system. It 
is possible that for certajn velocity, frequencies may become zero 


causing buckling instab ili'ly*' 

Per pinned— pinned and clamped— clamped rods one may regard 

the term Z ^ as a generalized centrifugal force; when this force 
9x^ 

overcomes the flexural restoring forces, the rod buckles. Por oanti- 
levered rods wi14i tapered free end buckling ii®^y occur when the lift 
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force generated at the inclined tapered free .end overcccies the restor- 
ing flexural forces. Although these effects are modified by the p 3 ?e— 
sence of viscous forcesj 13aey remain dominant so far buckling insta- 
bilities are concerned, Paidoussis [ 17 ] . 

3.4.2 FliDTlER IBS'PA'RTT.T'PTTCR . 

Oscillatory instabilities can be explained by conside 3 ?ing 
the energy balance. IShen the energy transfer from fluid to rod and 
vice-versa, in the course of one cycle of oscillation, exactly balance, 
the condition is said to have neutral instability, i.e., dynamic 
equilibrium. But when the energy transfer from fluid to IJae rod is 
33}ore than the energy transfer from rod to the fluid, the amplitude 
increases without limit causing flutter. In the opposite case osci- 
llations are damped. 

Paidoussis and Issid 27 have shown that conservative 
sjT'stema can only be subjected to coupled - mode flutter type of ins- 
tability, i.e«, frequencies of two modes are equal, in addition to 
buckling type of instability. But author finds that consenrvative sys- 
tem cetn also be subjected to uncoupled flutter, see Figs. 3.2, 3.3, 

3.4, 3.5 and 3.6. It is left for future workers to explain this. 

For detailed study of mechanism of instability one may see 
Refs.t 14, 17, 27 3 . 



CH/LP!I!SR VI 


C0I13IUSI0IIS 


In this chapter, the summary of the analysis made and the 
conclusions arrived at are illustrated as follows! 

1 • The finite element models developed for the new and old 

fCEniiula,tions are correct as the results obtained for the above 
mentioned formulations matdh quite well with liiose obtained 
earlier by classical methods for the cases of pinned-plnned rods 
and cantilever rods with tapered free end. 

However, for certain ranges of the flow velocities, ill- 
oonditicning is present in the system for the pimed-pinned rods 
and cantilever rods with tapered free end. Por these velocity 
ranges, the results obtaioed by finite element method and cla* 
ssical method may not be reliable « 

By increasing the number of finite elements the results do 
not converge for ill-conditioned velocity ranges for the cases of 
pinned— pinned rods and cantilever rods with tapered free end. It 
is observed that in general the Imaginary part of the comploc 
frequency converges while the real part shows little divergence* 
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The column, matrices { N. 2 m_„ V 


(e) 1 li 




.(e) ill 


ft \c;/ JLx I 

al— } 

,3,.(eP ,4 


.1e)V. h 


} and 


{ W^(EI — + lil — ) y found in Eqn* (2,14) » help 

3 3 1 x=L 

in Introducirg natural boundary conditions explicitly. These 
must be included in massj damping and stiffness matrices. These 


column matrices can not be neglected as suggested by Szabo and 
Iiee [ 243 , 


3. Por rods with pinned-pinned ends, coupled-mode flutter 
occurs for one panrticular velocity only, i.e., at critical velo- 
city for coupled-mode flutter, beyond which ordinary flutter 
occurs. IPhis contradicts the theory presented by Paidoussis and 
Issid [27] that apart from buckling, the conservative systems 
can only be subjected to coupled-mode flutter. It is left for 
the forthccsning researchers to explain. 

4. The results for clamped- clamped rod for old formulation 
are corrected. Early results obtained by classical methods , 
Paidoussis [ 14 1 , show the absaioe of buck l i n g in the first 
mode. V/hereas the results obtained by finite element method 
indicated that buckling will occur in tiae first mode and the 
critical velocity fcr buckling hsss been found to be very close 
to the first buckling velocity for conservative systems obtaingS 
by usirg Euler’s method, Paidoussis and IssidC 2? ] . 
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Sesults for clamped-c lamped rod for new formulation are 
plotted. The behavior of the clamped-clamped rod changes a lot 
by the inclusion of transverse component of tangential frictional 
force. Flutter instability does not occur in this case, only- 
buckling instability occurs in fiist and third mode. 

6. The finite element method vdta G-alerkin’s approach proves 

to be very versatile and flexible, because it gives a strong 
means to account explicitly for the natural boundary conditions 
and geometrical boundary conditions. 
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(a) Diagram of a clamped-f ree rod with tapered 
end in parallel flow 


Element numbers 



( b) Rod divided into m elements with m +1 nodes 


y 
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imension less complex frequency of a pinned-pinned red 
tion of the dimensionless u, for f3=0'1, 













/ I m (oo) 

The dimensionless complex frequency of a pinned-~pinned rod 
function of the dimensionless velocify u, for B “0*48. €c 1=0*2 5, 














6-55 EFFECT OF INTERNAL DAMPING 



the dimensionless complex frequency of a pinned -pinned rod a 
d function of the dimensionless flow velocity u , for B-O-l 
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APIEIDIX A 

BQUATIOIT OF SmiL lAIPEBAL MQglONS 

The system to be considered for the present study consists 
of a flexible cylindrical rod which has circular cross-section. The 
rod is immersed in a fluid which is flowing with a uniform velocity 
V, parallel to the neutral axis of position of rest of liie rod vdiich 
coincides with the x-axis of the co-ordinate system considered as 
shown in Pig. (2-a)» Cross-sectioiial area of the rod is Aj linear 
density (mass per ixnit length) is m^ » with flexural rigidity EIj 
and totsO. length of the rod is L, The density of the flovving fluid 
is P^. 

Per the formulation, follo\7ing idealizations are made; 

(i) The rod is assumed to represent a Bemoulli-Buler beam. 

(ii) The fluid flowing past the rod is considered to be incom- 
pressible. 

(iii) Small lateral motions of the rod about its positicn of rest 
are considered, during which the angle of incidence & and 
r— remain sufficiently small so that ; 

3 X 

(a) no separation occurs in cross-flow, 

(b) the fluid forces on each element of the rod may be 
assumed to be the same as those acting on a correspond- 
ing element of a long straight rod of the same 
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cross-sectional area and inclination. 

(iv) The fluid is supposed to be contained by boundaries suffi- 
ciently distant from the rod to have negligible influence 

on its motion. This assumption is made for the system which 
consists of a rod immersed in an iifinite fluid and is not 
valid for arod-shell system. 

(v) The virtual loass is equal to P^A for unconfined flow, where 

is the density of the material of the circular rod; and 

otherwise C PA, where C is coefficient of virtual mass, 
m r m 

Paidoussis I171 and Chen { 3] • 

(vi) The fluid flow is steady and uniform. 

(vii) The rod is supported at one or both of its extremities, and 
it is assumed that the supports are such that they generate, 
at most, only localized disturbances to the idealized flow 
conditions assumed here. 

A.1 PORCES TO BE OOESIIEEEI) : 

We now consider a small element 6x of the rod, as shovai 

in Pig. (2.a), where T is the tension, Q is the shear force, M is 

2 2 

the bending moment, m^( 9 y/9 t ) &x is the inertia force of the 

rod. P, and F, are the frictional forces per unit length in the 
IT jj 

transverse and longitudinal directions, respectively, and 

the rod gravity. Apart from these forces the rod is subjected to 

fluid pressures and hydrodynamic force which will be discussed in 
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detail along w/i'Ha the fricticmal forces in the next section. 

A. 1.1 PIIJID PR E SSITEE : 

let F and F be the resultants of the steady-state 
joc py 

pressure p acting on the outer surface of the element. These 

resultants can be found as follows : Consider the element 5x of 

Pig, (2-a), momentarily frozen and immersed in fluid on all sides 

as in Pig. (2-b); then there will be forces pA and pA. + » 

on the two forces of the element in addition wilii P^^ a^id ^py» 

The resultant of these additional forces and of P„ 5 x and 

px 

P 6x is equal to the buoyancy force. Assumiiog p = p(x) is a 

py 

linear function of x (which is reasonable to assume, as it applies 
to both a hydrostatic pressure distribution and to one modified by 
slcin friction pressure drop), it can be written (see Pig. (2-b)) 


EP = 0 

X 


P 6 X - (pA. + 6 x) + pA. 

px 9x 

= buoyancy force 

- A 6x 

3x 


ZP s 0 

y 


P 6x - (pA.+g^^) + H 6^) 


^ _ 


+ 

hence neglecting the higher order terms one obtains 
- p__- = ^ » 


ax ax 


(pA. *-^) ’ 

■“py “ax ax' 


(A-1) 
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Thus, from Eqn. (A-1) it is obvious that P is equivalent 

py 

to B tensile force 9 acting at th.e cross-section of the rodj of magni- 

g A 

tude pA* Also 5 if — = O9 1 = Oj hence for a rod of unifom 

o ^ 

cross-sectidi, the only pressure force acting on the element is equi- 
valent to a tensile force equal to P . 

py 

A, 1.2 HYDRO PIHMgC PQRCE : 

The inviscid hydrodynamic force P^ , v/hich is equal to the 
lateral force per unit length acting on the rod, and is given by the 
change in momentum of the lateral flow about the rod, Iiighthill [^IcT] , 
The components of fluid velocity in x- and y-directions are, 
respectively, 


V = Y cos 0 , 

X ’ 


Y = Y sin © > 

y 


(A-2) 


where 6 is the angle between the axis, of rod and the flow direction. 

Assuming small transverse motions of the rod in either direc- 
tion of its axis, 0 may be approximated by 

© = ( 1A) f ( 3 y/ 3 t) + Y ( 3y/ 3x) ] • (a-3) 

Therefore, fromEqns, (A-2) and (A-3) one obtains 


and 


Y = Y as cos 0=1, 

V = [ (3 y/ 3t) + Y (3y/3x)] . 

V 


(A-4) 


Momentum of the lateral flow of the fluid is 

°m ^ 8 y/ J » 



66 


hence ohe rate of change of momentum per unit length becomes 


(3 / 3t + V 3 /3x) (3y/3t + Y3y/3 x)l , (A-5) 


where m^ is -the virtual mass per unit leng'ih of the rod, which is 
given by 



(A-6) 


where D is the rod diameter, p „ is the fluid density# C is a ■. 

f m 

"virtual mass” coefficient vdaich is discussed by Chen [ 3] in detail. 

According to Idghihill [ 10 1 , the rate of change of momen- 

tum per unit length gives rise to an equal and opposite lateral force 
on the rod, i.e., 

^A = ^f ^ ^ V 9/3 X y. (A-7) 

A.1.3 PRIORI OUAl P QRCSS ; 

Friction forces acting on a rough cylinder set obliquely to 
a fluid stream are taken as proposed by Taylor [*25 J • For rough 
cylinders , he suggested the following expressions for normal and- 
tangential drag forces : 


®’lT = 1“ Pf D (C^ sin Q + Cjjp sin^ i;, (A-8) 

? ^f ^ °f 

viiere "q is the angle between the axis of cylinder and flow direc- 
tion and and C„ are drag coefficients associated with skin 

f RF 

friction a nd pressuare, respective]^. Hoemer's [8] expressions 
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of viscous fcsTces 8,lso reduce to Eqns* (A— 8) and (A— 9)j Paidoussis 
[17] . 

Por small oscillations, as mentioned previously, the angle 
of incidence Q can be approximated by Eqn. (A-3), 


e = (iA) t ( ay/a t) + v ( ay/Sx) 1 , 

hence Eqns. (A-8) and (A-9) reduce to 


and 


P = 
N 

P. = 




2 ^f 

1 D v2 


(a-io) 


A. 1,4 VISCOUS -DAIifflKG PORCE s 

The force representing the viscous damping effect can be 
expressed as, Paidoussis [jl7]] , 


P = p D C ^ 
D 2 f D 9 1 


Therefoore , Eqns. (A-10) may be vnritten as 


P - 
N 


P . D V c„ (^rf + V ^ ) + 

f f 9 1 9 X 2 


D 0. 


lx 

Da t 





(A-11) 


^•2 EQUATION OF ; 

The forces and moments acting on a small elemeht, 5x , of 
the rod, which has undergone a small lateral displacement y(x,t) have 
been determined and are shown in Pig', (2. a). The force balances in 
the X- and y-diroctions, and a moment balance, Paidoussis [17 I are 
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^ g + Ft + (F„ + 5’*) ^ - S’ 
3x r ° L 3J A 3x px 


= 0 


(A-12) 


3x IT 

9M 


- F„ - F. + P + F ^ ( T 9^) — m = Oj 

A py Igx 3x ^ 3x^ r » 


Q = - 


3 X " 3 X 


(El A) 3 (I 3^), 


3x 


'*3x^3t 


(A-13) 

(A-14) 


v/here the fact that y and its derivatives are small has been utili- 
zed, and it has been assumed that Euler - Bernoulli beam theory is 
adequate to describe the motions of the rod| moreover, a Kelvin-Voigt 
type of damping in the material of the rod has been assumed to apply* 

Substituting Eqna» (A-1) into Eqns. (A-12) and (A-13), one 

obtains 

^ (T + pA) + (m^ g -|^ A) + Ft. + (F^ + PjH = 


3x 


I- - F„ - F. + F. 


jj " A"3 X 


(A-1 5) 


3x 


l5+lr { (T + m 0. 


N ” "A "1 3x 3x ^ ^ 3X^ “r 


(A-16) 


Now, if y'^»0(e), then F^^ and F^ are of the same order , 
Paid ous sis I 17 1 , so that 

■*' ^N^ ( 9y/ 3^) “^0 

neglecting these terms in Eqn. (A-15) shfl- integrating it from x to 
I), upon combining it with Eqns* (A— 15), (A-1 6) and (A-14), one 


obtains 
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Ul - + El + P 

9 x^ at 


(T + p/l) + m g (L - x) 
i» r 


-|^ A (I - x) + J P^. dxl ^ + Pjj + g - 1^) 

X Jl ^ X 


^=0. 

St^ 


(A-17) 


The terms P. PI. and P,. hawe already been determined in 

il 5 JS ii 

sections A. 1.2 and A. 1.3* Now the terms (T + pA)^ - 

determined considering the case of an isolated rod in parallel flow* 

A. 2.1 C/.SE OP M ISOLfOIED BOD IN PARAUEl PIOI^ : 


Vl/hen the cylinder is isolated in a fairly large flow channel 
the pressure drop due to flow may be considered to be negligible, 


§JP._ p g . 
3x “ S * 


(A-18) 


Consideriiig the term. (T + pA)^^ , if ftis rod is free to slide 
axially (or is free completely) , at x = I, one obtains = •“ 

(since pA arises from the pressure cm the sides of the rod) , where 

p^ is the base pressures hence (T + p*l)j^ = Cpj^ - P^) 
p < p generally, this represents, more or less, "foim" drag, more 

accurately refe 3 ?red to as base drag, which may be expressed as 




(A..19) 
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If the total length of the rod is fixed between fixed supportss 

1 2 

at X = L there will be a compressive load - -r P D V C . (l/2) 

-K gj viiere the existence of the term g(li/2) depends 

on whether the length was fixed while the cylinder was horizontal 
or not. This fact can be explained as follows: 


Eqn, (A-I 5 ) can be rewritten after rearranging and neglect- 


ing the higher order term as 

A- + pA.) = - U g - ^ a) - . 


(A- 20 ) 


Integrating Eqn. (A-20) between the limits x to It one obtains 
tensile force distribution along the whole length of the rod as 

(I + pA)^ = (T + pA)^ + g •* A) (L - x) + - x) - 


(A-21) 


Now the stress - strain relation gives the strain in the rod 


due to tensile forces 


du- . 

dx ~ Fa ^^x’ 


(A- 22 ) 


where Uj^ is change in the length of the rod. 

Substituting Eqn, (a- 22) in Eqn. (A-21) and integrating the 
resultant expression between 0 to L, one obtains 

li II 

^0 I (T + + (m^ g - A) (1 - x) 


+ ( I - x) 3 dx 
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or 0 = ~ I (T + i + ® ** If - |-) + P^d-^ ~ f- ) 3 • 

.. (A-23) 

since fixed length requirement impli^ that uj, must be equal to zero. 
Solviiag Eqn. (A-23), one gets 
(T + eA)j. = - Pj^(V2) - (% S - g A) (V2) 

= -t Pf B Oj + g A] (1/2) . (A-24) 

Moreover j the re will also be a compressive load due to radial 
contraction, for a thin tubular rod this is equal to - 2V ( ^ ), 

And there may also be an externally imposed uniform axial tension 

T. 

Therefore for this case 

(T + pA)j_ = -[ -I D + m^ g - g A ] (l/2) 

+ (1 - 2 V ) ( ^ ) + I . (A-25) 

Combining Eqns. (A~19) and (A-25), and substituting into Eqn. (A-I?) 
along with Eqns. (A~18), (A-7) and (A-1l),C!ae obtains the equation 


of small lateral motions 
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+ El + C m ( 5T + ^ y 


.4 


8x 9t 


3x" 


-{ 6[¥ + (1 - 2V ) ( ^ ) 1 + 1 “ p^ 3) 7^ C^ + (m^ - mp g ] 


2 

[ (1 - ^6 ) L - Xl +1 (1 -6) \}‘H+1 Pf °f 

3x 

2 

( ^ , Y — pDCC l4-+(m~ Ei„) m ““p = 0, 

^ a+ + '' a-sr“ 2 D 3t ^ r ®3x r^u.2 


3-t 
(A-26) 


where 6=0 signifies that the downstream end is free to slide 
axially j 6 = 1 if the supports do not allow net axial tension* 
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ahekdix b 

E QUATION OB ffiTION (OID F O BMUIATIOir) AED I5JS PIKITE EIEi.ElIT 

ANALYSIS 


Neglecting gravity, pressurization effects, internal Ham p-i ng 
and the y-component of tangential frictional force , i.e., 

Eqn. (A~26) becomes 


4 


El + C (|r + V y - {6 ? 


X 


m f ^3t 




+ [ I D 1 [ (1 -6 /2) I - xl } 


3 X 

+ 5 If + v|i) + l 


+ m — ^ s= 0 

^ 9 t 


(b-1) 


Eqn. (B-1) represents the equation used by Paidoussisl 14 1 
which has been referred as "Old PormulatL on” in the present work. 

It does not take into account the frictional forces correctly* Com- 
pariiig this equation with Eqn. (A-26), one observes that in new for- 

afy 

mulation the term obtained by combining the coefficients of x 

and -If?' is P_ (x — “?• + (neglecting gravity) and in the old 

ox 9^ 3 V 

formulation -the above term becomes Ej.(x — ^ + 2^ ). In the new 

formulation integration by parts of the above termo is straight for- 
ward, see Eqn, (2,1 5), Integration by parts of the corresponding 
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term for old formulation vras not so straight forward. Ihis diffi- 
culty was overcome by splitting this term and then integrating by 
parts as follows: 


^ ^ = ^T. ( X ^ 1^) + • 

3 X 


■L-3^2 3.^ 


’Lax 


(B-2) 


Integration of all the terms of old formulation is given in 

3 V 

Chapter 2, except of the term • Applying finite element based 

Galerkin’s method, this term in Bqn, (b- 2) turns out to be after 
integrating by parts 

,(e) , h m. 


! "i *= 




0 


sf \ /"V . (B-3) 


Substituting Eqn. (2,25) into Eqn. (b- 3), one obtains for the element, 

U) 


K, Pj, y' 


V 


^4 V 


(e) 


h 


[k4] 


{y} j 


(B-4) 


where 


= F, 


:4] 

o 

H 

L 

. If J dK 

“-1/2 

h/10 

1/2 

-h/10 

-h/10 

0 

h/10 

-hVeo 

-1/2 

-h/10 

1/2 

h/lo 

h/10 

h^/60 

~h/10 

0 


(B-5) 
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and the boundary condition to be ccnsidered beccsnes for the Tdiole 


domain 


1/2 D y^®^ 


(B-6) 


+ 1/2 D y^®^ 


It should be noted that the values - 1/2 D 0 ^ and 1/2 D 0^ 

are summed up with ^ and 2 ^^ elements of the assembled 

stiffness matrix, respectively. 
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APBEHDIX 0 


BOUITOARY OOHDITIOIS FOR A CAM!IIEVt;?~ D HOD 'WITE TA'PRT^'R'n HREE BHD 


If "tile do’OTistrsain end is freej ‘teimina'ting in a "tapering endj 
"the cross-sectional area cf V\3aic]i varies smoothly from A to zero in 
a distance 1 ( << L), y and the lateral velocity v may be consi- 
dered constant for the length 1, Paidoussis [H] . This requirement 
allows the forces acting on the tapered end to be lumped and consi- 
dered in appropriate boundary c auditions, Hawthorne [ 7] . In the 
transverse direction equilibrium equation over the tapered end, say 
for L-l:x<L,is 


L 

J 


L-1 


3 Q 
Dy? 


dx 


f 

c 


I 'si 


j m (x) dx = 0 (c-l) 

L-r ^ 3t^ 


The parameter f , whi.ch is equal to unity for slender 

c ^ 

body, im'-iocid flow theory, has been introduced because the theore- 
tical lateral force at the free end may not be fully realised as a 
result of (a) the lateral flow not being fully two-dimensional, since 
tlie fluid has opportunity to pass around rather than over the tapered 
end [ Umlc, M.M., 1924 MCA Hep* Ho, 184 ], and (b) boundary-layer 
effects, Hawthorne [7] • Accordingly f^ will normally be less 
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than unity, Paidoussis [ 14] « 


Writing la^(x) = ^ A(x) and ia^(x) = ^ A(x) and, since 

y and v are constant over the tapered end, Eqn, (0-l) leads to 


I • - f 'V "V J ' 

I^^9x c f X 


2 L 


dx - ( p + f p ) — ^ / 


which yields 


where 


A (x) dx = 0 


+ PI — - + f ni„T(^'-^ + v|^)- 


(G-2) 


3x 3 1 


c f 3 t 3x 


(m + f Bi_) xe = 0 o 
^ st^ 


1 ^ 

7 ]■ A(x) dx . 


(C-3) 


(C-4) 


PurthcOTore , it may be assumed that there is no bending 
mcment at the free end, thus 


El + P I " 

3x ■ 3 X 3 t 


(0-5) 
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APIEKDIX D 


THB GOMPUTER gROSBAp.^ 


The computer programme was prepared for assembling "bhe 
system matrices, substitutiiig the boundary conditions, inverting 
the assembled stiffness matriic, forming flie dynamical matrix and 
finally for finding out the complex frequencies of the system. 

The computer programme was wri't't®^ ^ Portran 17 and run 
on an IBM 7044 computer with apprcxim£i:ely 32K words of core sto_agei 
The prograime is quits gemral ara built-in subroutines ttllSV for 
matrix; iOTSrslon and ZBIDIA for findtas out ocnplex frequencies were 

called directly, 

Fortran listing of the program iis Appendrs 


D,1 deSCBTPTIOH OB ERO&BAtJ : 

The MA-HT programme 3 ?eads and prints various input d 
output results, other subroutines calculate the matrices, asseubln 
fhem and after applying boundary oonditicns foim the dynamical matr 
Ihs oompler elgenyalucs axw calculated by calling built-in subroutine 

a INDIA . 


S ubroutine Na me 
STOT 


Description 

Called by MA.IN, calculates the 
various elemental stiffness matrices, 



Subroutine 


JIIASS 


DAMP 


ASMBL 


BNDHY 


MTim 


Name 


Descriptlcai 

total stiffness matrix, assembles 
tliem, Introduces boundary conditions 
by calling BilDEY and calculates 
inverse by callii^ MATIW, 

Callad by MIN, calculates the ele- 
mental mass matrix, assembles and 
introduces boundary conditions by 
calling ASMBIi and BKDEY subroutines* 
Called by MIH, calculates the ele- 
mental damping matrices, sums them 
up, assembles them and introduces 
boundary conditions by calling ASIBL 
and BKDBY subroutines in turn. 

Called by tlASS and BALIP, assembles 
the mass and damping matrices in 
turn. 

Called by STIFF, LIASS and BMP, 
introduces boundary conditions to the 
assembled stiffness, mass and damping 
matrices and separates the required 
matrix. 

Called by MU, foms laae dynamical 
matrix’. 
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Descripticaa 

Called by STH'I', computes matrix 
axiverse « 

Called by MAIN, determines eigenvalues 
and eigenvectors of a general real 
matrix. The binary version of this 
programme is stored on 1301 disk file, 
which is named Z INDIA (Nor further 
details refer to communication of 
ACM journal, Vol 2, 1968, pp. 820)# 

LIST OP ERINCIB IE VARIABIES 

PCETRAN Implementat ion D efinitim 

Program Sym bol 

IIAIN mOGRAM 
AL 
AMI 
AMP 
AMP 
AREA 
AV 

BETA 

CBAR, CDASH 
CM 


length of the rod. 

Moment of inertia of the rod* 
Mass of fluid per unit length. 
Mass of rod per unit le33gth 
Cross-sectional area of the rod 
Actual velocity of the fluid 
Dimensionless mass of fluid 
Priction drag coefficients. 
Coefficient of virtual mass* 


Subroutine Name 
MATINV 

EIGENP 



FORTRAN Impleaentatim 
Program Symbol 


Definition 
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PROGRAM 
DO, DI 

E 

evin(i) 

bvrn(i) 

FAC 

FD 

I 

NMIN 

NMAX 

NSP 

OMEGA1, 01/BGA2, 0]ffiGA3 

ROW 

RO»P 

VNON 

SVNON 

StTBROUTIMB STIFF 
AXP 


Outer and inner diameters of ttie rod, 
respectively. 

Modulus of elasticity. 

Imaginary part of dimensionless 
eigenvalue. 

Real part of dimensionless eigenvalue. 
Dividing factor viiich gives dimen- 
sionless eigenvalues. 

Tangential drag force. 

Number of elements. 

Minimum number of elements. 

Maximum number of elements. 

Total number of velocity points'. 

Lowest three exact natural frequencies. 
Density of fluid. 

Density of rod. 

Dimensionless velocity. 

Dimensionless velocity data. 

Axial distance from upstresaa end to 
node p, see Fig, (l-* o). 
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11 

P 

GSB 

GSM 

S 

51 

52 

53 

54 
S6 
I 
Z 

SIBROUCDIHE MASS 
M. 

GMB 


Definition. 

Length of one element. 

Largest element of GSM 
Factor multiplied to elemental 
stiffness matrix XHDHD (4j4) 

Size of the assembled matrix =r 2 x (l+l) 
Size of the separated matrix 
Factor multiplied to elemental beam 
stiffness matrix 
Assembled stiffness matrix (N,i) 
Separated stiffness matrix 
Total elemental stiffness matrix (4>4)« 
Elemental beam stiffness matrix (4 $4). 
Elemental stiffness matrix MFVSQ ( 454 )'* 
Elemental stiffness matrix DY/DX (4,4)« 
.5 lemental stiffness matrix XKDIL (4j4)* 
SI + S3 - S4 

Factor multiplied to matrix S3* 

Factor multiplied to matrix S2 

Elemental mass matrix (4j4)* 

Assembled mass matrix (NjN). 
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Program Symbol 

SUBROUT H® JMSS 
GMM 
H 
P 

SUBROUT 11® DAMP 

Aim 

C 

01 

C2 

03 

GOD 

COM 

N 


P 

PBi.Kr 


SUBROUTli® jlSkl^L 


AA 


ASM 


N 


Definition 


Separated mass matrix (lTl,ll), 

Size of assembled mass matrixs:2x(l+1) 
Factor multiplied to matrix AM. 

Factor multiplied to matrix 02 
Elemental damping matrix AMPAY ( 4 , 4 ), 
Elemental damping matrix DY/DT ( 4 , 4 ), 
Elemental internal damping matrix (4*4) * 
Total elemental damping matrix ( 4 , 4 ), 
Assembled damping matrix (n,1T), 

Separated damping matrix (iTIjll), 

Size of assembled damping matrix 
= 2 X (l + 1). 

Factor multiplied to matrix 0, 

Factor multiplied to matrix 01. 

Elemental matrix (4j4)» 

Assembled matrix (I 5 N). 

Total number of degree of freedom 
of tbe syststt (size of ASM)=2 x(i+1). 
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Prograni Symbol 

Definition 

SUBROUTIMl BKDRY 


ASM 

Assembled matrix (HjH)* 

BSM 

Separated matrix (Hi jH1 ) 

U 

Dimeiasion of ASM = 2 x (l+l). 

N1 

Dimensicaa of BSM. 

SUBROUTINE MATMIY 


D 

Unitary matrix I y (h1,H1) 

Dl 

[k"”* m] , (hi, HI) 

D2 

[k“'* cj 5 (hi, HI) 

DMM 

Dynamical matrix (HB, HB). 

HI 

Dimension of j^Kj , |^m] , [^Cl , 


[l] «d[oj 

HB 

Dimensim of dynamical matrix = 2xHl 

SUBROUTINE MATINV 


A 

Coefficient matrix of siae Hi to be 


inversed . 

B 

Constant vector of size Hi 

M 

Humber of constant vectors 


(if Ms 0, only inveise is computed)* 

mmm 

Value of tbe determinant 

SUBHO\iTli® EIUBWP 


A 

Matrix of size H1 x Hi for whicb. 


eigenvalues and eigenvectors aw 
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Prograxa Symbol 
SI3BROUTIHE EIGBNP 
EVI (l) 

EYE (l) 

EM 


T 


YBCI (J,l) 
VBCR (J,l) 


Imeginary part of eigenvalue 
Real part of eigenrvalue 
Dimension of dynamical matrix giren 
in the main programme 
Number of binary digits in the 
mantissa of a single precision 
floating point number = 27,0 
Imaginary part of eigenvector 
Real part of eigenvector 
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C 
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■ CWP IS IFF DFNSITY OF THE ROD. 

•/ur-i IS TFE DIMENSIONLESS FLOW VELOCITY. 

' VNCS IS THP DIMENSIONLESS FLOW VELOCITY DATA. 
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/ ■■{ ",>:C[,-AP = *,5: £.e,5X,*CC4SH=*, „16.S,5X,*EETA^ 

r = ' '.]"=■ fly F)/( 1 , -e'^'TA) 


;i 6 «P) 




( 




♦■• c 


)/( :-. 1'I*(CCI**2) 5 


5X, *FA C = 


f ( 4fM )/{ { ^yFS«CM + AyF)*( AL.**4 5 5 ) 

: y] AN F, FAC, Af £fl 

vi = -S"X,Fie.e,5'-:,*AMP=*,t.l6.8 ,5)',*AyF=*,F16, 

J i ‘ " V .lA-.F) 

■ ’‘rr'-L^ 7 * r ^ c 

p c c 

' ■ ^ i' ■ = .• a , C * r; ^ r, 

• ■ " ' f , ■ f' ^ ' fl ’ , I M='(: *, £A3 

. - 5 ^,‘''Cf'?:GAI=*,Ei6,P,5X,*jMFGA2 = *,?16,B,5X,*!'MrGA3=*,ei6, 85 

' fli'^ V.CN/ .4, 3,'", 2*''),3,0,4,0, F.u, 6, C , 7, 0 , 8. (J, 9. 0 , 1 D. 0 , 1 1 . D, 12 ,'f 

" 2 P~ 1 *- 

; " 5 ’,.' Sr »i« X' 2 i: 7 li ^ 3 ^ 3 ^ 3 ^ afc^e sjs ^aS* 4 ^! 4 ' 3 ^ 5 §'i 35 fi# 5 Ss 4 «aif# 56 «#sj 5 #jfic 3 ^:^ 3 ^ 3 |«+#^-' 

■„i I fv r: , 

, y;. M' :/ '"i ■ 

■" :■! i-t yiN,5MAx 

Mi :.'' jjj = i,fiSP 
v'‘,T.:‘rTOV'.jfi!.{ JJJ 5 

^ V= { GM^'T ( (F'^'A^I ) /AMF)*VNO' i) /AL 
f Li( FCWF*f.C*( AV**2)*CBAR)/2, 

^'RnT77,I,VKrN,AViFL 

CALL OTIFF 
C a L L y A S 0 
CALL DAMP 

call ^'AT^'FY 

n I Ft- ^'S I ON CVR( 30 5,£VI( 30 5,VECR{30,30) ,VECI ( 30, 30) , I NDIC ( 30 5 
OIRf'NSinN FVRN{ 20),EVIN'( 30) 

T*27,0 

'!1={2*( I + l)-2) 

’! = ?*(;! 

NH = 3 « 

CALL PIGFNP(K,NM,DHM,T,eVR,eVI,VECR, VeCIfINDIC) 

" ' ' * 3 ' 'f''' ' ' ''L 'I. '''' N''' '' ’' ' ' ■ ■' 

n£AOF=rlFVP( L)**2)+(EVI{L )**2) } 

fVR{ L ) = EVRf L l/DENOW 
CVI(L)--F:VI(L>/D6NCH 
rVPM(L)=EVR( L)/FAC 
''V JN{L) = EVI< L)/FAC 

COM INUE 




c 

. , ..... „. ...v- 

:!i j'“. 3‘s .1:’' 

'4-* ■•■•A# 4U-r %3U 

.,.,lt J-'l* .ifi-i'* W^\ 


. ^ ^ r •' r , r 

iT- 

TIFF 

/•' 


A- :’.i *^1 

.Ti*. .#■>(». *R|-5 

1'* 

'. • ' .i. ■ i ■ > 

1. ": 

.^lAL CC 


:. i 'T r 

"'• 1 ... . ' . ■ -0 

"Fl; 

L -NC-TF 

'C' 

''': T'*' 1 

r 7^, 

■: LA ROE 


'■ ; ,,.,4 „ 

T f" 

.H rL'Y 


, , ' |... 

■'■' 1 

. ML - TiL 

c. 

. T r -r ^ 

- :.[ 

‘VENTyiL 

f’ 

. i E ■" F 

■•:■ .-,L 

- lvc.‘..T *1 

r,' 

... T 

::. ;■'■ | 

f-’C* T^L 

r. 

, , ., ,« 

H ' 

c 3 r ./ ^ L F- 

r 

. T 


■" .3 • A J C 


- 

/ " f 

1 :,) 


" ■ '■' ' r ■ / ' 

., ^ - 

.''( 2., ,2" 


■■ .t-y 

i 

'<'■ F , A V 


'. /:: 

^ .1 / A 

0 I » E 


■ ./!: 

’ 5 / A 

F, AMP, F 


./■_ 

■ S / 

.1 


'- : ,L t 

1 T / T 

, JJJ 


;■ 

1 1 / c 

. { 1 ,.‘,10 


r ■ (v,..r . / 

:• ' / :' 

'{ 10,1C 


■ '5 / . 

~ ^ / n 

.',DI 


r , /f. 

"‘I li 

•=( 1 i,i.r: 


^ 1" / 5 

"T/S 

A( Iv-i 30 


' I y -.■ N 5 I il 

'0 S^ 

{4,a5,g 


'■■| I y » r, I "' 

' p 

(2'', 15 


: 35* ^ ^ :<^35s 3§! 4« sfe 5|« ::^ ^‘c^ 'it 3 §c :i§i -ifr ajt '4* ^ ^ # 4« --> 


•'’n!Di:->/?.T" OF THE iLE.Mt <T. 
,F C''ie EL-iKEiT, 

,T EL:.ME'“T DF matrix GSf^« 
..>TAL STIFFflESS MATRIX. 
50AM STIFFNESS MATRIX. 
STIFFNESS f'M'-nX MFVSO, 


ITIFFOESS HATRIX DY/CX, 

;ti ff u- ss mate ix XNnoo, 

STIFFNESS MATRIX. 
STIFF’M'SS MATRF y. 


'IR-MsSIiOO SMAtA),GSB( 2 C, 2 G) ,GM 8 ( 20 , 20 ) ,GCB(?.r), 2 C) 

AIR - * SI-'- 

-)IR NSI 5 K‘ AA(l*',l!!),ASM{ 20 , 20 ),BSM( 20 , 2 C) ■ 

■■ AL ■I.R*'',ALtAMI 

Q j*i :$y/ y- 1. * ?;< * >!■ =<« * ^ M y' 5> * Xi sj # sJi* X" ####*** ‘4# Xta? ♦ ♦***!{:*^ 4t #!>!##!!'*♦ Xs ** ♦♦ 

r»AT 4 CN, CSLTA/l.R, i.O/ 

DATA Cl-:/. . / 

r.ARflrl >. ■ 

'.= {CAKA>.‘J~x<Ay n/( ( fiL*X■3)=X(AHP-A^^F) ) 

*T 15A, C-AYAtG 

fopj^ATf /, 5X, *GAMA=*,F16.8,5X.r«G=>x,E16..8) : 

F y« yi S)i ys yr SJi *:{: y: 

nc,aMA=D« I . ■ ' 

TFfli: = (CGA(^A4.t=!'AMI)/{AL**2) 

pr; INTX 5 S ,CGARA,TeAR 

l '»5 FCPMAT(! 5 X,=!.CCA^^A=’t',::l 6 , 8 , 1 CX,= 5 ‘T 8 AR==X .£ 16.8 ) 


156 

+ 

c 


HATA PTA,Mj/D.0,0.3/ 
PAEAP,= (ETA*E*AHI )/{AL’»X<2) 


PRIIvTl56,FTA.NUf PAEAR , 

Fr:RMAT{5X,*F«A = #,F16.8,5X,«Nl’=*,El6.8*5X,*PABAR==f=,Ei6.8) 

ys yt :$ ± :!^ yt y? * # jJ; yt 4 ys 3}! aj! * yt ^ # 1# 4 * 4 x. # ♦ * # * # 35!!ie 5(5 ★ ♦ St: 4. Xt Jj. ♦ 34« * * * ♦ ¥ # .«I ♦ * # >Jt ijr + 3)C * :#(♦ # 3(! ♦ sjt 

ffLPT = AL/FLCAT'(I J 

4'. Jif yi * )!i 4 * sj. .«. * « ♦ * * 5jf sJ. 4 4. ♦ * « # * * * * * * ♦ # * 4t 45 sSs Xt afc # * # ♦ Jft ♦ J(t # sje * Jfe * # ♦ # « * 4 3«i ♦ * ♦ ♦ # Jji ♦ 


p - ( f:«WI ) /( ELHT’*>^2 ) 



S{1,1) = 12.«P , ; 

S ( 1 ,2) = 6,>^ELPT«P 
S ( 1, 3l=-5( if 1) 
r> l!,4) = S(l,21 
S (2,2 )*4,M.(ELyT**2l«p 

S.{2, 2|=-SC 1,2) 

S (2,^)s52.*|f l,MT45f2)*P 

DO i ^.''orrr-. ;,'v" . ‘ • 


' I,' t f. , 

r;'-.. M/OO',;;. 


Millii 



r'»4 

?3 

HI 


121 


r I 


4 I 

V / 


5 F )=i'C 


» * f -' 4| J 

i,-! I 


- ) = 

) 

A ) . 

' )-(■'! " ')=»!<? ) 

:)=“5^( :,2 ) 

F )=- ( 

)=^Al 1 ,1) 


#^'4N 


'■ ' ( , 

' 

■. r ‘y 

{ 


; K= 1 , 4 


,0' = :’ ,4 

■ ■, ( y'' 

, K ) = s 

•*■ r D 

1 «t 

« » L ^ 

' ■■ -j 




"•*,(■( C 4 


■■, f 

r.“ 

■ ;i ■'- 


( / 

u 


f IK , JK ) 

.JJJ*GF,3) GO TO aO 


1 *iX 1 A : . <; « P ) 


1 ( ,. 

3(1 
1 ( 1 

1 (’ 
1(1 
i C’ 


1 ,1 


. F>''X 

, f L'^TAX ) /ITi, 

, 1 ) = > . .F ’I' X 
, A ) =- S 1 ( 1 1 2 ) 

,1 )=-Sl( 1,2J 

,F) -'• '■ 


> 


1 ( 2 t M = S 2 ( ! , 2 ) 


1 {' ,^ )=- (X=^ (ELf'T**:!) ) /6e« 
i (3,1)=--C.5’5<X 
•>! (3,7)=-Sl( 1,2 ) 


GO TO 81 
MATRIX DY/OX 


IS’*') 


)=-Sl( 1 . _ . 

l(3,3.)=-Sl( 1,1) 

F3 t2,A)a=:Sl (1,2) 

FI (4,2 ) = S1( 1 ,2) 

Fi {4,2 3=’-Sl( 2,4) 

Fl{A,3)=-Si( 1,2) 

FI {A,Al=f'«7 

IF(I.Nt.2,CP.JJJ.GE.3) 

7 PINT 1 <4 

FORM AT { / 1 5X , =^51 1 FFN8 S S riA i m ut/ua 
»PIMT?3, { (S1(IK, JK),JK=1,4),IK=1,4) 

FrjPMAT(/^4X,4E20,8) 

COM I NUN 

•■1 = 2* (I + l) 

DO 5 IL=I,4 
DC 5 JL=1,4 
r.6ill, JL)=O.D 

S6(U, JL) = S6(IL, JL) + {SiIL,JL) + SllIL,4L)-S4f IL, JD) 

CCNTIf^UE 
nn 121 iK=a,N 

DO 121 JK=l,K ^ , 

r,S«3( IK, JK)»0,0 
A.XP'''3»y- . ^ 


hiiimm 

AXP»AXP4i| 



)*C 




, 5* 

5(1 


■ ' - , ? :* *( ‘V 



:\:h y{X 


*'3';'tr 

j4i«, ! 



» ' 


f ]: 


, 1 ) 

} -■■-:?( -.r 2) 

J ~ i f 

T^')=li2 flK, Jl') 

_ t 3) GC T ■' 79 

■f/,;x,''‘5TiFFiN.,SS MAT!?!)?' KFVSC iS>!'} 
' T i ( n ■ ( IK , JK ) , JK=1 ,4),IK=1,A) 


; I'' 


i » 


{ 


J‘'=l , 


J<) = G^ { IK, JK)* S 2 {IK, JK) 


'! : ~ 


,C'- .JJJ 


•*2) GO TO 82 


'■A 


■ t/, 5 X,ATCTAL STIFFM 2 S 5 MATRIX IS*) 
'-'i , ( (33 ( IK, JK ) , JK = 1,A) , IK=1,4) 

■( / , AX, A 2 '.8) 


9 


37 

8.3 


13 


106 

38 

EA 


30 


! J=-I I 
I I 

J J3 = ^ J+'3 

I?''’ IK = I! ,113 
')•. LII JK = JJ,JJ3 
L = IK~ Il4l 
L J=JK™JJ+1 

OSH IK , JK ) = GSB( I K, JK) + S3( L,LJ) 

OOMIfJJi- ■ . 

'ALL BNGRY{ GSBtGSM) 

# 3 ( 1 : :j|!3{c4E^**9K* ******** ******** ********** 

■■I = i2=^( I+l )-2 ) 

FACtr'’=ABS( C-SM{ 1,1)) 

.7C 9 IK = ],h 
OC 9 JK = l,fi 

I F (AES (GSM IK, JK )).GT* FACTOR) FACTOR -A8S CGSH (IK , JK ) ) 

r.CMINUF: 

IF{ I,NF.2,(;R.JJJ.GE.3) go to 83 

'’RI'.TITjFACTCr 

FORMATf 5X,*FACTDP=*,E16.8) 

COh'TlNUf- ■ . ■ 

nC I 3 IK = 1,K 

r)Oi'’jK=i,) ■ ■ 

nSMflK, JK ) = GSI«( IK, JK) /FACTOR 

COM I hue 

|F{I,KE.2.CP«JJJ,GE.3) GO TO, 84 

PRIhT IJ6 

F0FMAT(/,5X, ♦ASSEMBLED STIFFKESS/FACTOR IS*) 

PRINT38, ( (GSM( IK,JK),JK=l,N),IK*l,N) 

rCPVAT|/,5x,eEl5-8) 

COMIKUE 

CALLING SUePOUTINE HATINV, 

no :3c;ii=i,N 

# 1 1 -0»p 


107 

39. 

' : '85 ' 




8( _ 

CALL PATIK' 

PRINT 10 
FOPMAtf/. 
PPIHT39»f 
,FCRMAT( 
iXONIlNUg 


DpSERM) 

GO to 85 

f y ' -I Om"' ^ ' 


M I 




UUUUiJU 


I 


L - ,JJ.J 

;= T r, ’ y • 


;) GO TO £6 


, ::vy OF GS.f^/FACTCP; i £« 5 

, f ( r- M KI J ) ,KJ = 1,; ) 

( / , : X , 0 , = ) 


r 

+ 


1.09 

2 

8 ? 


5 


-r ^ ^ 5 ^ 3 ^; jjs t? Jr 5 t« 5 ^ S^ti^S ^ '# 3 ^ 3 ^ 5 $< 1 ' 3 ^ ^ V ^ 

' h t s '^' ' ' ' ' 

Jj 4 ' 4 "J. ^ sji :$! 3 §s 4 ' 4 : 3 ^ 3 ^;^ sgs ^ tf,s 3 ^ sjc j}s # 3 ^ 3 |s ^ # sir 4 : 4 f t 3 ^' * 


]••: t;- •-L’l’^'iiaL «ASO 'MATRIX. 

r. T r , r c C(.'D|I ^^ 5 $ p ft T R I >{ , 

’ *’': "'F-' O'PE-AT'C MASS MATRIX. 

•; ,' •-, :'::>*!'■!> 5 i" ij- :< Si 3 f. sjt ,V : 5 ^ ^ $ :jr 3 {t !{i :};<! jjc^: ?[; j 4 : 3 fe!j 5 Sf; sjt sj: ;(; 1 ^ si! 

f' - ' /•„'.'/F f ;- -, 1 . ) 

f.«r { 2 ' ir *' 5 

-' '“ / l]/AfT,AV 

/.-: I '/AFtAPF^FL 
/i. :F/AL 

/-'I5/!:CWP,FCWF,APEA 
/,.l 7 /r,JJJ 
T"N GPF (2 »F 0 ) 

ICO AA{ ic,ic' ),ASP{ 2C,2'. ),&SM( 20,20} 


?fs 5^ 4! 3^ ❖ ♦ ^ # *4t # ^ afs s!? sji 1,; -A* 


" V r > 

• N ^ f * 

” ;’ \ K f 
- ^■- \, A r ’ 

■A i ^ : ^ ’ 

a: 

' \iil 

L^'Tr/iL/FLC^TJ I ) 

a: 3^' 7/ A's 4 «,r 5?^ 4 3l« ^ 4 3^ 4 444 44 4 3 ;c4 44 44 44444444444 44444444444444444444 44443|4444<4 

n = ( { APF + APF)=>ELMT)/a 20.0 

'•• ( 1 , n = 1 9 6 . « F 

•■1 { 1 , ? ) = 22 .*ELMT=i'P 

■4 { 1, 3 y = 5'4 • 4 F . .'■■ ■, 

•‘■i ( ] , A y=-i 3 .*r;LPT*p 
M ( 2 ,?)= 4 , 4 ={frLPT’!==> 2 )=i'p 
M{?, 3 )=-P(l, 4 ) 

M( 2 ,A)=- 1 .«(ELPT=i>* 2 )*P 

M{ 3 , 3 }=M( 3 , 1 ) 

M ( 3 , 4 )=-M( 1 , 2 ) 

“( 4 , 4 )=P( 2 , 2 ) 

no i 11 = 1,4 

nc 1 JL=l ,4 
H ( JL,II ) = F( II,JL ) 

I F{ I.NF. 2 .CP.Jj}J.GE. 3 ) GO TO 87 

PRINT 1(19 

FORMAT!/, 5 X,«WASS MATRIX IS>»=) 

PRINT 2 , HM! II, JL ), JL = 1 , 4 ),! 1 = 1 , A) 

F 0 RMaTl 2 X,AF 20 . 8 } 

continue 

CALL ASMBLfM.GMB) 

CALL BNCRY(GHB,GMMy 

RETURN 

ONC 




REAL .•EtMTi.AI,. AMFiAV ^ , . .. 

C'' IS ^TK '6 ' ELEKENTIi; DAWPlfl© HATRIX , 

SA|RIX« ' 

Sp'' ' ■ '■ 

I#.' >. 




n ' 

'-in 

d 3 


C 


lU 

31 

89 


-, ' r ' / . 


- -f / 


■ : 1. - /:. 

' / i P , W F 

- ■ ’■•; / : 

■ 7/ I , JJJ 

■' / : 

/CiA'NCDA, 

‘- 

'’3/C!(l 

' ■' /,. 

■A/‘':(i <,’C) 

' C / ^ 

/ c c , 1 1 


vie (2 >j?.’i) 


’ ' AMI 


* j. 

{ I ) 


f 


1 i 

■ » 


]=' 

) = 

) = 

) = 
)- 
) = 
)r, 

1 = 

\ =. 


i 

"S ?*' 

'm\/i . 

R 

"C,*£LMT-':'C 


5 

* 5-: 


^ : 5 e 5 ^ 9 ^i|e 4 cfjs 5 <*s s|e:^«isjs :fr 5 |« j|e 


3 §r »!' 


p 


(.[ M ##;’)* B 

? ' . e 

... ..,*i: 

) = ( ..LVT >>=^2 ) *2 
)=- 6 ,*RLKT >!‘0 
) = '..»• 

.%%2.CR. JJJ,GE®3) 
ll'- 

C/^MPING MATRIX IS!^} 
JK= 1 , 4 )» IK= 1 , 4 } 


I f 


C ( A - ^ 


G ( A 1 A 
I F f I . 

f'PIXT 


GC TO 88 


rf’PKAJI / j fX^^FIRST CAMPING 


OHMAT 1 / , f I 

..r,INT3y., ( (£.( IKtJKK. . 

FCPMATf /,AX,AF20.8} 
r/'ihT INUfi 

•Ie aH ♦ i # 9ff # s^:§t ^ 4 5^3^s§! ^ sj« ^tafs 3^aSta$t3|3$e3|s3®6s|c3|e%55c#3ie35c3|c5|s 3|k3^*:^2|!2|j5^3^aS«3|!-»i}« 

PI^KI = { {RCWF*(CC/2. )*AV*CBAR )+CDASHJ /420.0 
Cl f l ,1)=156.*ELMT’*'PINKI 
C3 (I I ?) = 22.’^ (ELMT’!'*;:)’5‘PINKI 
C1{1,3) = 5A.*ELMT*PI?4KI 
CU1,4 3=-12,#(ELMT**2}t=PINKI 
C 1 { 2 , 2 ) =4 ( PLMT>1'^3 ) =^'PINK I 

Cl [2,3)=-Cl{ ]»4) 

Cl f2,4}=~2.*{kMT**3)=«'PINKI 
Cl (3,3l=ClCl|l) 

Cl (3,4) =-C IE 1,2) 

C: |4|4)=CH2,h 

^ «!f # itt 4i ^ 9{> :i^ 3(; :$ 3(1 3)i 4: 4 :$i 3^ # ji> # 3» 4 ^ 9t( 4c « aft # 4 ^ 4 a|c # 4^ 4c :4c * # 4: 3^ * 4= 

DC 4 1K=1,4 

DC 4 JK=1,4 , . 

Cl (Jf<,IK) = ClClK, JK) 

t p / f .♦■ir . 9 .re - I I j -CP 


■!..» ■, ■Wf- :f .* i:^ f — V * .VI 

IF(I,.NE..2,CR.*JJJ..GE.3) GO TO 69.. 

^PINT 111 

FORMATE /,5X,*SeCOND DAMPING MATRIX IS*) ' , 

PPlNT3i,( (CHIK,JK),JK®1,4),IK*1,4) . 

FCPMATE /,4X,4e20,8) ' 

CONTINUE . .. . . . .'..'''' ......... . ... ... ........ . . .-..' . . .. ..’i:,-.,a „'.■....., 

;44443jc444444 44!4444c«44 443»,3|c443(:3)t4i44[4:4)44;4c](:«4i444444;44c3(t4tiic4i4;4c3|;44>***4f****i*! 

^lUhT-r’"’ -•' 


1350 - 


T ,, 






oon o ^ ooo 



' jj: 5*; “4f iiV :‘f ^ :*c ^fc 3^ sjc ^ :4s sjt ije 3 ^ 3^.::^'4 sjr '3^ 3 ^ %.3^ ^ 3 ^ 5?c •;§£ 4c 3§: :^! .:|£'5|c 4c# 3fii3|c 

'■ tr cijt r-:;; iSvfiL(A^,/iSP) 

!i{ij},;!4j(s|i*-J:J4#!t!ji:J:4r3^'}i43|jii^ 4 4-.4I4: 4< * 4' =4= #’♦’}'#♦ #’5' 

',A n Th'" lu-nsm^l '/atrix, 

’ sv Ij THc aSGEMBLED matrix. 

^ r.CNMCI'.'/Hl 7/1 , JJJ 

■> IM':;hr>IC"i AA(10,10),ASM{2C,2C),BSM(20,20I , . 

'! = ?*{ I + i) " ^ - 

TC I 11=1 ,V' ■ ■ , :■■ • . A : 

'AfJ 1 JL=1,N ■ ■; - , 

1 -tSf^l II, JL) = ''t«'> 

nc 2 K=!,i ■ ■..:■•■■. ;■ 

A I I = 2’!'K--i 

iL=n ■ ',■ -'A' ^ 

I 13=11 + 3 

■ Jt3=JL + 3 ■. 

OQ .? IK=IUn3 • 

■irj 2 JK = JL, JLl . . 

L = IK-n + l \ ' . 

LJ = JK”.JL+1 

2 ASH IK, JK ) = ASM( IK, JK) + AA{L,LJ) 

I F(I.f,r:.2.CR .JJJ.Ge.3} GO 70 79 

PRINT i''’J3 

1<A3' FdRMATt'/,5X,=J'ASSEMBLE0 MATRIX AFTER B. C. IS#I 

PRINTIO^, ( {ASM|IK,4K),JK»l,H),lK»lfN) ■' ' 

J.:'?4 FGRMAT{/,5X,6Ei5.H 
79 CChTINUE . 

PcTUPN . ... 

FNC ^ 











9 


'r i 11 = 1 


. .JL=%K 
( T I» JL 

n-, (II^'jlUc! f II , JL ) + (”GSf^(II,K)*G^'M(K, Jl H 

T iN’jf' ' ■ ■ ■ ' ' ■ 

FnKwp,9''»^TRIX C2. ' 

'U, F IL=1 ' 

2 JL=1,K 
•?nL,JL) = C,f 

HL, JL3 = C2(TL, JL 

■ r I. T T f ^ 1 1 


) + (-GSM{ IL,L J’S'GCMJ L, JL)} 
MATRIX .Df'M. , 


-'1 

'’={2’=>( I4l )-? ) 

!fc = 2=^^ ■ 

'■I!;' 1 IK-1,NP 
"iC 3 JK = I, Ne 

{ IK, JK)=C.fl 
l=^+l 

DO 4 IK=L,^e 
D.P 4 JK = 1,K- 
I I = IK»N 
JL = JK 

OMK{ IK, JK) = DNM{IK,JK) + 01{ II, JL) 
DO *1 IK = f. ,KB 
DO 5 JK=L,KB 

JLaJK«^| 

D^f'dK, JK)«O^MtIK,JK) + D2I IL, JLI 
DO 6 IK=<f1,^ 

' _JK*L,WB ■ 

•OR. 4 JJ. 66. 31 60 TO 91 
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